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ABSTRACT 

Context. The Sun shows abundance anomalies relative to most solar twins. If the abundance peculiarities are due to the formation of 
inner rocky planets, that would mean that only a small fraction of solar type stars may host terrestrial planets. 

Aims. In this work we study HIP 56948, the best solar twin known to date, to determine with an unparalleled precision how similar 
is to the Sun in its physical properties, chemical composition and planet architecture. We explore whether the abundances anomalies 
may be due to pollution from stellar ejecta or to terrestrial planet formation. 

Methods. We perfom a differential abundance analysis (both in LTE and NLTE) using high resolution (R ~ 100,000) high S/N (600- 
650) Keck HIRES spectra of the Sun (as reflected from the asteroid Ceres) and HIP 56948. We use precise radial velocity data from 
the McDonald and Keck observatories to search for planets around this star. 

Results. We achieve a precision of o" < 0.003 dex for several elements. Including errors in stellar parameters the total uncertainty 
is as low as cr 0.005 dex (1%), which is unprecedented in elemental abundance studies. The similarities between HIP 56948 
and the Sun are astonishing. HIP 56948 is only 17±7 K hotter than the Sun, and log g, [Fe/H] and microturbulence velocity are 
only +0.02±0.02 dex, +0.02±0.01 dex and +0.01+0.01 km s"' higher than solar, respectively. Our precise stellar parameters and 
a differential isochrone analysis shows that HIP 56948 has a mass of 1.02±0.02Mo and that it is ~1 Gyr younger than the Sun, as 
constrained by isochrones, chromospheric activity, Li and rotation. Both stars show a chemical abundance pattern that differs from 
most solar twins, but the refractory elements (those with condensation temperature Tcond ^ lOOOK) are slightly (~0.01 dex) more 
depleted in the Sun than in HIP 56948. The trend with Tcond in differential abundances (twins - HIP56948) can be reproduced very 
well by adding ~ 3 Mg, of a mix of Earth and meteoritic material, to the convection zone of HIP 56948. The element-to-element 
scatter of the Earth/meteoritic mix for the case of hypothetical rocky planets around HIP 56948 is only 0.0047 dex. From our radial 
velocity monitoring we find no indications of giant planets interior to or within the habitable zone of HIP 56948. 
Conclusions. We conclude that HIP 56948 is an excellent candidate to host a planetary system like our own, including the possible 
presence of inner terrestrial planets. Its striking similarity to the Sun and its mature age makes HIP 56948 a prime target in the quest 
for other Earths and SETI endeavors. 

Key words. Sun: abundances - stars: fundamental parameters — stars: abundances - planetary systems 



1. Introduction 

In recent years there has been an important number of 
studies related to solar twins, stars which are spectroscop- 
ica lly almost identical to the Sun. The reader is refereed 
to ICavrel de Strobell (Il996h for a review of the early his- 
tory regarding the search for solar twins, that culminated 
with the identificatio n of 18 Sco as the closest ever so- 
lar twin (Porto de Mello & da Silval 119971: ISoubiran & Triaudi 



iden- 
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l2004l) . More recently, new solar twins have been 
tified iMelendez et al., 2006; Melendez & Ramirez, 12007 : 
Takedaetal., 2007; Pasauini et al., 2008; Taked a & Taiitsu . 
2009; M elendez et al.. ,2009: .Ramirez et al.. ,2009h and HIP 
56948 has dem oted 18 Sco as the star tha t most closely re- 
sembles the Sun dMelendez & RamfrezLl2007l : lTakeda & TaiitsuL 
2002). 

Solar twi ns are useful to calibrate the zero-point of the 
temperature ("Casagran de et all 1201 Ot 'Melendez et al.', '2010B|; 
Ramirez et al., 2012) and metallicity (Melendez et al., 2010B|; 
Icasagrande et al.L 120 111) scal es, to better char acterize the in- 
teriors of stars like the Sun (iBazot et al.l l201 li) . and to iden- 
tify the transport mechanisms that cause Li depletion in 



the Sun (do Nascimento et al"?, '2009; Melendez et al., "^2010^ 
Bauman n et al., 2010; Denissenkov, 2010; Castro et al., 20lj|; 
ILi et al.l 120121) . But most importantly, solar twins are the per- 
fect targets to look for small chemical abundance anomalies that 
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may have been tirinotic ed in previous works (lGustafsson[|2008t 
iGustafsson e"tani2010l) . 

The first spectroscopic study designed to exploit the advan- 
tages of differential abundance analysis between solar twins and 
the Sun, showed that our Sun has a peculiar chemical abun- 
dance pattern, suggested to arise from accretion of material de- 
pleted in refractory el ement s due to the formation of terres- 
trial planets (.Melendez et al.L |2009|) . Further studies have con- 
firmed, with different degrees of accurac y, that the Sun i ndeed 
has an anomalous surface composition ('Ramirez et al.', '20091 
I20T O; Gonzalez et al., 2010; Gonzalez Hernandez et al., 201ofl. 
We show in the appendix A, that the reality of these abun- 
dance anomalies is well establis hed. Besid es the important im- 
plications for planet formation (IChamber s, 2010) and to ex- 
plain abundance anomalies in Jupiter (iNordl und. 2009), the solar 
abundance peculiarities may be re levant for modeling early stel- 
lar evolution ( IBaraffe & ChabrieH l2010h and for solv ing the so- 
lar modelling crisis when using low solar abundances (iNordlundl 
l2009tlGuzik & MussackLl20Tol) . 

Even before the deficiency of refractory elements in the so- 
lar convection zone was discovered, Castro et al. (2007) investi- 
gated the effect of accretion of metal-poor material onto the Sun 
as a way to help solve the solar modelling problem whe n a low 
oxygen abundance (jA llende Prieto et al., 2001; Asplun d"et al.L 
12004 ': 'Melendez & Asplundj, 120081) is adopted, and they found 
that indeed accretion pr ovides some improvemen t, but the prob- 
lem is not fully solved. iGuzik & MussaclJ (1201 Ol) demonstrated 
an improvement in the comparison between stellar models and 
helioseismic data when accretion of low-Z mateiial is taken into 
account (their Fig. 9), albeit a full resolution of the discrep- 
ancy is not found. More de t ailed modelling has been recently 
presented bv lSereneUi et al.l (l201 lb . who use up-to-date nuclear 
cross-sections and include accretion of metal-poor and metal- 
rich material, considering a range of accreted mass and differ- 
ent timings for accretion. They conclude that there is somewhat 
better agreement with helioseismology for differentiated accre- 
tion, but not complete agreement. Overall, models with metal- 
poor accretion improve the agreement with the helium abun- 
dance inferred from helioseismology, while metal-iich accretion 
improves both the depth of the convection zone and the sound 
speed profile, with exception of a model with late accretion of 
0.015 M0 of metal-poor material, which improves the agreement 
with the sound speed profile. 

A detailed test of the terrestrial pla net formation hypothe- 
sis was performed by IChamber"sl (1201 Ol) . who used the compo- 
sition of about two dozen chemical elements in the Eaith and 
CM chondrites (representative of the asteroid belt). Interestingly, 
IChambersI (120101) showed that Eaith material alone can not fully 
explain the peculiar solar pattern, but that a mix of Earth and me- 
teoritic material gives an excellent fit for more than 20 chemical 
elements. Thus, the peculiarities in the Sun could be a signature 
of both the formation of teiTestrial planets and of the asteroid 
belt. 

An interesting alternative interpretation of the solar abun- 
dance anomali es, from an analysis of the solar twin M67-1164 
(lOnehag et all 120^1 1). is that the chemical peculiarities may re- 
flect that the Sun was bo m in a massive open cluster like M67. 
Using various arguments, lAdamsl (120101) concludes that the birth 



' I Gonzalez Hernandez et all ( l201(]h have contested the pla net sig- 
natur e scenario, but further scrutiny of their work by Rami rez et al.l 
demonstrated that in fact the results of Gonzalez Hernandez et 
al. are fully consiste nt with the works of iMelendez et al.l ( l2009h and 
iRami'rez et aP (l2009l) . 



environment of the solar system might be a moderately large 
cluster with lO'^-lO"^ members. Nevertheless, only one solar twin 
in M67 has been analyzed to d ate for high precision chemical 
abundances (lOnehag etal.Ll20rH) . so, more observations are ur- 
gently needed to verify if indeed all solar twins in M67 have the 
same solar abu ndance pattern . Also, notice that based on a dy- 
namical study, iPichardo et al.l (l2012h have shown that the Sun 
could not have been born in M67. 



IGustafsson et al.l (I2OIOI) warned about potential systematic 
effects in chemical abundances due to differen t fines of sights 
when the Sun is compared to the solar twins. iKiselman et al.l 
(2011) have recently studied the line-of-sight effect using fiigfi 
resolution observations at the solar equator and at latitude 45°. 
Seven key chemical elements in a broad range of co ndensation 
temperature were analyzed bv lKiselman et al.|(l201 ih . who show 
that there is no difference in the abundances obtained at differ- 
ent latitudes for both volatile (to within 0.005 dex) and refrac- 
tory (to within 0.002 dex) elements. Thus, it is very unfikely 
tfiat tfie abundance anomalies seen in the Sun ( Melendez et all 
2009; Ramirez et ak, 2009) can be attributed to line-of-sigfit in- 
clination effects. 

In appendix B, we show that the Sun's chemical peculiari- 
ties also do not arise due to the particular reflection properties 
of the asteroids employed in the analyses, as expected given that 
the relative reflectance of asteroids show mostly smooth changes 

19961: 



1995 



Binzel et al.l 



over hundred s of A (e.g. iXu eta 

IBus & BinzelL 120021: itazzaro et al L 120041: DeMeo et all l2009l) 
about 3 orders of magnitude wider than the narrow stellar lines 
used in abundance analyses. 

Only some 15% of solar type stars appear chem i cally s imi- 
lar to the Sun (IMelendez et al.L 120091: iRamfrez et all |2009|) and 
therefore we regard the Sun to have an anomalous chemical 
abundance when compared to other solar twins. Assuming that 
the solar abundance anomalies are due to terrestrial planet for- 
mation, then perhaps only these 15% of solar type stars that are 
chemically similar to the Sun may host rocky planets. This is a 
lower limit to the amount of rocky planets formed around other 
Suns, as part of those planets may fall into their host stars, alter- 
ing t hus the original abu ndance signature. The Kepler mission 
(e.g. iBorucki et al.L 12010 ) will give us the first estimate of the 
frequency of Earth-sized planets in the habitable zones of solar 
type dwarfs. 

It is important now to find, through a detailed chemical 
abundance analysis, stars which are chemically identical to the 
Sun and which may therefore potentially host other Earths. HIP 
5694^ is the perfect candidate for identifying subtle chemical 
anomalies, as this star has been found to be the most s imilar 
to the Sun in stellar p arameters ([Melendez & RamfrezL 120071: 
iTakeda & Ta]itsuLl2009l) . The first high precision (cr ~0.03 dex) 
detailed abundance analysis of this star showed that HIP 56948 
may be one of the r are stars with a solar abundance pattern 
(iRamfrez et al.L l2009h . In the present work, we perform a much 
more refined study (cr ~ 0.005 dex) of HIP 56948, to assess its 
similarity to the Sun and to which extent it may host a planetary 
system like ours. 



a.k.a. Intipa Awac/ian http://tiemeylab.blogs.nytimes.com/?s=intipa| 
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2. Observations 

2. 1 . Keck HIRES spectra for high precision abundance 
analysis 

HIP 56948 and the Sun (reflected light from the Ceres aster- 
oid) were observed with HIRES ( Vogt et al., 1994) at the Keck 
1 telescope on May 19, 2009. Exactly the same setup was used 
for both HIP 56948 and Ceres, and the asteroid was observed 
immediately after HIP 56948. A total exposure time of 800 s 
was used for both HIP 56948 and the Sun, consisting of mul- 
tiple observations co-added, in order to avoid non-linearity. We 
stress here that for highly accurate work an asteroid should be 
used instead of the daytime skylight, as there are important vari- 
ations in the skylight spectrum with respect to the solar spectrum 
(iGrav et al.L |2000). Furthermore, asteroids are essentially point 
sources for typical observing conditions (seeing > 0.5 arcsec), 
thus the observation and data reduction for both stars and the 
asteroid are performed in the same way. 

A resolving power of R x 10^ was achieved using the E4 
0.4"-wide slit, accepting some light loss (seeing was ~ 0.7 arc- 
sec) in order to achieve the necessary spectral resolution. For 
HIP 56948 the signal-to-noise level measured in continuum re- 
gions is about 600 per pixel at 6000 A, while it is somewhat 
better (S/N = 650) for Ceres, for which the predicted magni- 
tude at the time of the observation was V - 8.3C B i.e., some- 
what brighter t han H IP 56948 (V = 8.671±0.004. iOlsenLll993t 
lRamirezetani2012h . 

The spectral orders were extracted using makeeQ. Further 
data reductions (Doppler correction, continuum normalization, 
and combining spectra) were performed with IRAF. A compari- 
son of the reduced spectra of HIP 56948 and the Sun is shown in 
Fig-H As can be seen, the spectra are nearly indistinguishable. 



2.2. McDonald and Keck radial velocities for planet search 

Soon after HIP 56948 was identified as the best solar 
twin jMel endez & Ramirez, 2007), the McDonald Observatory 
plane t search program (e.g., lEndl et al.L l2006t iRobertson et al.L 
l20T2h began to monitor this star The observatio ns have been 
carried out with the Tull Coude Spectrograph (Tull et aO Il995h 
at the 2.7m Harlan J. Smith Telescope. HIP 56948 has so far 
been observed nine times (from May 2007 to March 2012) and 
the scatter of the observations is 5.2 m s"' (discarding one out- 
lier), which is consistent with the 5.1 m s"' median error bar 
The radial velocity data is presented in Table [Tland Fig.|2] 

In February 201 1 we started observing HIP 56948 for planets 
using HIRES at Keck during time allocated to a NASA key sci- 
ence program to support the CoRoT mission. We have acquired 
30 datapoints with HIRES up to February 2012. The radial ve- 
locity data is presented in Table|2]and Fig.|2] 

For both instruments we use a temperature-controlled io- 
dine cell for wav elength calibration. We use the Austral code 
(lEndl et al.L 120001) for the computation of precise differential ra- 
dial velocities. 

The scatter of the HIRES observations is 4.4 m s ', which 
is higher than our 2.8 m s"^ median error bar The difference 
amounts to 3.4 m s"' and can be explained by typical jitter val- 



' http://ssd.jpl.nasa.gov/horizons.cgi 

'* MAKEE was developed by T. A. Barlow specifically 
for reduction of Keck HIRES data. It is freely available at 
http://www2.keck.hawaii.edu/realpublic/inst/hires/data_reduction.html 

' The HIRES data for HIP 56948 were obtained at times when the 
CoRoT field was unobservable. 
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Fig. 1. Comparison of HIP 56948 (red circles) and the Sun 
(solid line) in diff'erent spectral regions. The quality of the 
Keck/HIRES spectra is very high for both HIP 56948 (S/N ~ 
600) and the Sun (S/N ~ 650). It is hard to distinguish any dif- 
ference between both stars. 



ues as measured in other stars. For example. IWrightl (1 20051) finds 
a median stellar jitter of cr^j^ - 4.4 m s"' for inactive stars with 
(B-V) > 0.60. Interestingly, the McDonald observations seem 
to be less affected by stellar jitter, due to the much longer expo- 
sure times required when compared to the Keck observations (~3 
min), damping the radial velocity variations due to stellar jitter. 
Indeed, one of the observational strategies to reduce stellar noise 
in short time scales (mainly stellar oscillations) is to take rela- 
tively long exposures (e.g. a long 10-min exposure, or the aver- 
age of several sh ort exposures), thus imp roving the precision of 
the observations (iDumusque et al.L I2OI lb . We verified this with 
our Keck data taken during 4 nights in January 2012. Each night 
we took three consecutive short (~3 min) exposures, so that the 
average exposure for a given night has a much lower contribu- 
tion from stellar noise (lower panel of Fig. |2}. The night-to-night 
scatter is only 3.3 m s"', in reasonable agreement with our ob- 
servational error bar of 2.8 m s"' . The Keck measurements taken 
in February 201 1 and in January-February 2012 clearly rule out 
hot Jupiters, and we can even eliminate the presence of planets 
with masses as low as Neptune in the inner 0.04 AU (lower panel 

of Fig. la. 

Considering that velocity semi-amplitudes about two-three 
times the typical measurement precision (in our case ~4 m s"') 
can be detected with confidence, we can rule out the presence 
of a stellar companion and of nearby giant planets. As shown in 
Fig.|3] where the 2-cr (dashed line) and 3-cr (solid line) sensitiv- 
ity of our observations is shown (computed as the planetary mass 
that would introduce detectable radial velocity variations), there 
is no indication for a Jupiter-mass planet in the teiTestrial planet 
region (<3 A.U.), and even a less massive giant planet such as 
Saturn may be iTiled out (at the 2-cr level) inside 1 A.U. So, the 
inner region around HIP 56948 seems free from giant planets. 
Examples of some circular orbits for a Jupiter-mass planet at 
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Fig. 2. Upper panel: precise radial velocities obtained with the 
Tull Coude Spectrograph at the 2.7m Harlan J. Smith Telescope 
(HJST) of the McDonald Observatory (filled circles) and with 
the HIRES spectrograph at the 10m Keck telescope (stars). We 
show some circular orbits for a Jupiter-mass planet at 1 AU (long 
dashed line), 1.5 AU (dotted line) and 3 AU (solid line). Lower 
panel: average of three short (3-min) consecutive observations 
taken with HIRES/Keck during four nights in January 2012. 
Some circular orbits due to hypothetical Neptune-mass (solid 
line), Saturn-mass (dotted line) and Jupiter-mass (dashed line) 
planets at 0.04 AU, are shown for comparison. No giant planets 
have been detected so far in the inner regions (<3 AU) around 
HIP 56948. 



1 AU (Earth's distance from the Sun), 1.5 AU (Mars' distance 
from the Sun) and 3 AU, are shown in the upper panel of Fig. |2] 
We have verified for a range of planetary masses, distances, ran- 
dom observing times and including error bars in the measure- 
ments, that no giant planet is present, although notice that since 
our data are sparsely sampled, some giant planets (especially on 
orbits with e > 0) may have escaped detection. Once more radial 
velocity data are obtained in the coming years we intend to per- 
form detailed simulations to see which kind of planets we could 
have missed. With the existing dataset, we see no evidence of 
giant planets in the inner region around HIP 56948. 

Providing we maintain a radial velocity precision of about 
4-5 m s ' for our observations, in a decade or so we should be 
able to detect (or rule out) the presence of a Jupiter twin, i.e., a 
Jupiter-mass planet orbiting at 5 A.U. from HIP 56948. 

We will discuss the implications of the constraints we have 
deduced which limit the presence of massive planets in the inner 
regions around HIP 56948 further in Sect. 4, in combination with 
our findings of the chemical similarities between HIP 56948 and 
the Sun. 



3. Abundance analysis 

The abundance analysis is based on the Keck HIRES 
spectra, following our previous differential work on so- 



Fig. 3. Estimated sensitivity of our observations as a function 
of distance (upper panel) and orbital period (lower panel). The 
2-cr and 3-cr sensitivities are shown by dashed and solid lines, 
respectively. The red dotted lines represent planets with 1 Mjup 
and 1 Msat- Our radial velocity data discards planets as massive 
as Jupiter in the terrestrial planet region. 



lai- twins (iMelende z et al.'. '2006'; 'Melend ez & RamfrezL 120071: 
Melende z et al.L 2009; Ramirez et al., 2009), where the solar and 
stellar spectra are measured in exactly the same way. 

An initial set of equivalent width (EW) measurements was 
obtained by fitting gaussian p rofiles with the automatic rou- 
tine ARES (" Sousa et ail 12007'). We computed the relative dif- 
ference in equivalent width between HIP 56948 and the Sun, 
SWr = {EW* - EW^)IEW^, and fines with 6Wr deviating from 
the median < 6Wr > for a given species, were measured by hand 
both in HIP 56948 and the Sun. All weak fines (EW < 10 mA) 
and the lines of species with only a few lines available, were 
also measured by hand. About 20% of the EW measurements 
needed to be checked according to our initial empirical analysis, 
although note that this procedure only reveals the most obvious 
outliers. The deviating automatic ARES measurements could be 
due to any of these causes: /) the number of components found 
by ARES in the local fitting window is not exactly the same in 
the spectrum of HIP 56948 and the Sun, //) the 2nd-order poly- 
nomial used to fit the local continuum could be somewhat diff'er- 
ent in both stars, contamination by telluric lines do not fall 
exactly on the same place in both stars. Indeed, most faulty auto- 
matic measurements occur either in the blue, where the spectrum 
is more crowded, or in the red, where telluric contamination is 
higher. 

An initial model atmosphere analysis is performed using the 
improved EW measurements and again we look for outliers from 
the mean abundance for each species, and they are checked by 
hand. In most cases the revised measurements result in a reduced 
scatter. The improved EW is mainly due to more consistent man- 
ual measurements, which are performed using exactly the same 
continuum for the Sun and the twin, and also the same part of 
the profile, or the same treatment of blends. This is not always 
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the case with the automatic ARES measurements The final anal- 
ysis is performed using the set of revised equivalent widths, ex- 
cept for lithium which is analyzed using spectral synthesis. The 
adopted EW for HIP 56948 and the Sun are given in Table |3] 

We use ID Kurucz overshooting model atmospheres 

(ICastelli etall [T997I). as well a s MAFAGS-OS (Grupp, 200l 
and MARCS JGustafsson et all I2OO 8) ID LTE models. The 



models have diff'erent mixing-length approaches, as explained in 
detail in the above references. The differences in stellar param- 
eters, between HIP 56948 and the Sun, are small, therefore es- 
sentially the same results (within ~0.001 dex) are obtained with 
either Kurucz or MAFAGS-OS models. 

The analysis has been performed both in LTE and NLTE. 
For the LTE ca lculations we used the 2002 version of MOOG 
(ISnedenL IT973I) . while the NLTE calculations are described in 
Sect. 3.3. 

3.1. Stellar parameters 

The stellar parameters adopted for the Sun are Teff = 5777 K and 
log g - 4.44 (e.g. Cox, 2000). With T^ff and log g set, the micro- 
turbulence velocity (v,) is found by requiring no dependence of 
A(Fej3 with reduced equivalent width EW^ (- EW / A) for Fel 
lines. We found vf = 0.99 km s"' both for Kurucz and MAFAGS 
models. The above set of parameters (5777 K, 4.44 dex, 0.99 km 
s ') yielded the zero-point solar abundances for each line /, Af 
for a given model atmosphere. 

For HIP 56948, an initial set of stellar parameters was found 
in LTE using Kurucz models. The relative spectroscopic equi- 
librium was achieved using differential abundances (5A, for each 
line i. 



6Ai = A* - Af. 



(1) 



Thus, the effective temperature is found by imposing the rel- 
ative excitation equilibrium of (5A, for Fel lines: 



d(6A^'')/dO(,,,) = 0, 



(2) 



while the surface gravity (log g) is obtained using the relative 
ionization equilibrium. Usually, this is done using Fel and Fell, 
but in our case we verified that within the error bars the ion- 
ization balance is fulfilled simultaneously for Fe, Ti and Cr. 
Therefore, we use the mean relative ionization equilibrium be- 
tween Fel and Fell, Til and Till, and CrI and CrII: 



^Fell-Fel ^< ^^Fell > _ ^ ^^Fel ^ 
^Till-Til ^^Till ^ _ ^ ^^Til ^ 
^Crll-Cd ^< ^^Crll > _ < ^^Crl ^ 

A"-' = (3aP'="-F'=' + 2A™-'^'^ + A'=''"-'='"')/6 = 0, (3) 

with the weights arbitrarily chosen, but reflecting increasingly 
larger errors for Fel/Fell, Til/Till and CrI/Crll. 

The microturbulence velocity v, was obtained when the dif- 
ferential abundances 6A^'^^ show no dependence with reduced 
equivalent width EW,-: 



d(6A^n/d(EWr) = 0. 



(4) 



The spectroscopic solution is found when the three condi- 
tions above (eqs. 2-4) are satisfied simultaneously, and when the 



A(X) = log(yVx/A'H) + 12 



metallicity obtained from the iron lines is the same as that of 
the input model atmosphere. Notice from the equations above 
that our work is strictly differential, i.e., we do not enforce abso- 
lute spectroscopic equilibrium, which may be difficult to achieve 
(both in LTE and NLTE) eve n in the Sun dMashonkina et all 
1201 iHBergemann et alll20T2l) . 

The initial LTE solution with Kurucz models showed that in- 
deed HIP 56948 is extremely similar to the Sun, with differences 
(HIP 56948 - Sun) in reif/log§/[Fe/H]/v, of only 17 K / 0.02 dex 
/ 0.02 dex / 0.01 km s"' . We also tried the MAFAGS-OS models 
and the same stellar parameters were obtained, with a negligi- 
ble difference of +0.001 dex in [Fe/H] from Fel and Fell lines 
with respect to the Kurucz models, and with a difference in the 
spectroscopic T^ff of only 0. 1 K. 

Given the very similar stellar parameters between HIP 56948 
and the Sun, we do not anticipate NLTE considerations to pro- 
duce significant changes to the stellar parameters. Indeed, the 
NLTE Fe abundances indicate the same T^ff within 0.4 K, with 
the LTE temperature being slightly cooler. In Fig.|4]we show the 
individual iron abundances as a function of excitation potential 
of the Fel lines, both in LTE and NLTE. The NLTE abundances 
result in ever so slightly lower line-to-line scatter (cr - 0.009 
dex, s.e. - 0.001 dex) than in LTE. 

Unfortunately the trigonometric Hipparcos parallax for HIP 
56948 is not known with enough precision to provide better con- 
straints than our spectroscopic value. The Hipparcos value is log 
g = 4.37±0.07 (TablelH, which agrees within 1-cr with our re- 
sult from spectroscopy (log g - 4.46±0.02). The above error 
in the trigonometric log g is due to both the uncertainty in the 
Hipparcos parallax and typical errors in photometric tempera- 
tures (e.g. iMelendez et all I2010B). Adopting instead our more 
precise stellar parameters we obtain a Hipparcos-based gravity 
of log g - 4.37+0.05, which agrees with our spectroscopic grav- 
ity within 1.3-cr. 

Notice that the Hipparcos parallax (15.68+0.6 7 mas , ac- 
cording to the new data reduction by van Leeuwen, I2007h im- 
plies a distance of 63.8+2.7 pc for HIP56948. We can obtain 
an independent estimate of this distance assuming that the abso- 
lute magnitude of H IP56948 is identical to solar; i.e. My =4.81 
(iBessell etallfTm) . Since V=8.671+0.004 and the eiTor in the 
absolute magnitude from the uncertainty in our stellar parame- 
ters is 0.055 mag (see eq. 8), we derive a distance of 59.2+1.5 
pc. Thus, there is agreement within 1 . 1 cr for both distance esti- 
mates. 

The adopted stellar parameters, errors, and comparison with 
other estimates, are given in Table |4] Overall there is a good 
agreement (within the error bars) with other independent es- 
timates. The errors depend on the quality of the spectra of 
both HIP 56948 and the Sun. Since the S/N is very high 
(>600) we can put stringent constraints on the stellar param- 
eters. Nevertheless, we are also limited by the degeneracy be- 
tween Teff, log g, [Fe/H] and v,, which increases the errors. As 
described in the appendix C, for a fixed log g (or a small range of 
log g values), we could determine T^s to within 0.8 K, while if 
Teff is keep fixed (or within a small range), log g could be deter- 
mined to within 0.006 dex. Regarding the microturbulence, for 
a fixed Tgff and log g, v, could be determined to within 0.0004 
km s Due to the degeneracies between stellar parameters and 
to the observational uncertainties, actually the errors are consid- 
erably larger (Table 0, 7 K, 0.02 dex and 0.01 km s"' for Teff, 
log g and V,, respectively. A detailed description of the determi- 
nation of stellar parameters and the uniqueness of our solution, 
is given in the appendix C. 
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= 0.0206 (tj= 0.0094, s.e, = 0.0014) 



solar spectrum. Molecular lines of CN ( [Melendez & Barbu^ , 
1999'; Mandell et al., '2004) and C2 ('Melendez & Coheni. i2007 : 
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Fig. 4. Iron abundances versus excitation potential of Fel lines 
in LTE (red stars) and NLTE (blue circles). The blue solid and 
the red dashed lines show the fit in NLTE and LTE, respectively. 
The line-to-line scatter in NLTE is only cr = 0.009 dex and the 
standard eiTor is 0.001 dex. 



As a further check of the effective temperature, we have tried 
to compare synthetic Ho- profiles, computed using MAFAGS- 
OS models, to the observed Ha profiles in HIP 56948 and the 
Sun. Unfortunately Ho- falls on an order too close to the edge 
of the chip, making it difficult to normalize that region properly. 
Despite the above problem, our tests using Ha indicate that in- 
deed HIP 56948 is somewhat hotter (~ H-20+20 K) than the Sun, 
in agreement with what is found from iron lines. 

3.2. LTE abundances 

The adopted atomic data is presented in Table[3] Whenever pos- 
sible we use laboratory oscillator strengt hs, or theoretical gf- 
values normalized to laboratory data (e.g.. iFuhr & Wiesel 120061: 
[Melendez & Barbuv, 2009). However, the input ^/-values are 
not critical, as they cancel out in the line-by-line differential 
abundances M,. The interaction constants Ce were co mputed 
from the bro adening cross-sections calcul ated by Barklei ii et al.l 
(I2OO O) and Barklem & Aspelund- Johansson (2005), using the 
transformation given in Melendez & Barbu y ( 2009). If broaden- 
ing cross-sections were not available, we multiply the classical 
Unsold constant by 2.8. The adopted Ce values are given in Table 

m 

The mean < 5Aj > and the standard deviation are computed 
for each atomic species, so that we readily identify any outliers. 
The suspicious measurements are checked by hand, both in the 
Sun and HIP 56948, using insofar as possible the same measure- 
ment criteria, i.e., the same continuum regions are adopted and 
exactly the same part of the line profile is used for the gaussian 
fit. 

Equivalent widths were used to obtain abundances for 
all elements except for lithium, which was analyzed us- 
ing spectrum synthesis, as in our previous work o n so- 
lar twins (Melendez et^ l2006t lMel6ndez & RamirezL l2007t 
iBaumann et al.^ ,2010) . The line list used for spectral synthesis 
is presented in Table |5] The data for t he Li doublet was take n 
from the laboratory d ata presented in Andersen et al. (1984). 
Although [Smith et al.1 (1998) and Hobbs et al. (1999) reported 
new ^/-values based on theoretical calculations, the difference 
with our adopted laboratory values is only Other atomic 

lines near the Li feature wer e taken from Mandell et al. (2004) 
and the Kuruc^[ and VALD dKupka et all 2000) databases, ad- 
justing in some cases their ^/-values to better reproduce the 



Melendez & Asplund, 2008) were also included in the spectral 
synthesis. 

The LTE abundances were computed using Kurucz models 
and checked using MAFAGS-OS models. The agreement is ex- 
cellent, with a mean difference (MAFAGS - Kurucz) of only 
-0.0010 dex in the differential abundances, and with a element- 
to-element scatter in the differential abundances of only 0.00075 
dex (0.17%), which is quite remarkable and shows the weak de- 
pendence of our strict differential analysis to the adopted model 
atmosphere. [^ 

The LTE differential abundances are provided in Table 
[6] but note that the adopted abundances are those based on 
NLTE (when available). The observational errors, which depend 
mainly of the quality (S/N) of the spectra of both HIP 56948 
and the Sun, are given also in Table [6] The observational error 
is adopted as the standard error (- cr/ V«) when more than three 
lines of a given species are available. Otherwise we assumed a 
minimum value of cr^i^ = 0.009 dex (s.e. = 0.005 dex for 3 lines; 
s.e. = 0.006 dex for 2 lines), which is the typical line-to-line scat- 
ter for species with more than 10 lines available. When only 2 
or 3 lines were available we adopted the maximum value of the 
observed cr and o-„^in, i.e., maxicr, 0.009 dex). If only one line 
was available, we estimated the eiTor by performing a number of 
measurements with different assumptions for continuum place- 
ment (within the noise of the spectra) and profile fitting. 

In addition to the observational errors, we also present the 
errors due to uncertainties in the stellar parameters in Table [6] 
where the total error is also given. As shown in appendix B, our 
very small observational errors (~0.005 dex) are plausible. 

3.3. NLTE abundances 

For most of the elements we have been able to account for de- 
partures from local thermodynamic equilibrium (NLTE). For 
Li, C, O and Na we have employed MARCS model atmo- 
spheres (Gustafsson et al., 2008) an d the model atoms described 
in Lind et al. (2009), FabbianetaD ([2006, 2009) and lLind efsH 
(201 lA) . respectively, using up-to-date radiative and collisional 
data. We note that for Li and Na, quantum mechanical estimates 
of the cross-sections for collisions with both electrons and hy- 
drogen are available. For O we adopt a scaling factor Sh of 0.85 
to the classical formula of iDrawinl (I1968L 11969 ) for excitation 
and ionization due to inelastic H collisions as emp irically deter- 
mined from the solar center-to-limb variation by iPereira et al.l 
(2009). For C we adopt Sh =0.1 in the absence of similar 
empirical evidence. For these elements the NLTE calculations 
were performed with the ID statistical equilibrium code MULTI 
(ICarlssoni [19861) . 

In addition, we performed NLTE calculations for Mg, Al, Ti, 
Cr, Mn, Fe, Co, and Ba, usin g the revise d version of the DETAIL 
code (originally published in lGiddin"gsl 1198 1) and the SIU code 
(J. Reetz, un published). We used MAFAGS-OS ID LTE model 
atmospheres (Gruppi [20041) provided by F. Grupp (piivate com- 
munication). The differences in abundances required to equalize 



http://kurucz.harvard.edu/ 



The MARCS mod el and the mean atm ospheric structure of a 
3D model atmosphere dAsplund et al.L |2009|) . are much closer to the 
MAFAGS-OS model than to the Kurucz overshooting model, thus the 
effects of using either the MARCS or the 3D model would be even 
smaller than for the comparison between MAFAGS-OS and Kurucz 
overshooting models. 
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NLTE and LTE equivalent widths are defined as NLTE correc- 
tions. 

For Mg and Al, the model atoms were kindly provided by 
T. Gehren; those models were previously used in the spectro- 
scopi c analysis of solar-type stars given in iGehren et al.l (|2004 
| 2006l). Atomic models f or Cr, M n, Co, and Ti, were tak e n from 
iBeraemann & Cescutti (2010), Bergemann & GehrenI (l2008h . 
iBergemann et al. (2010), and Bergemann (2011), respectively. 
To compute NLTE corrections for the lines of Fe I/II, and Ba II, 
we constructed the model atoms from the laboratory and theo- 
retical data given in NIS10andKuruc2E] databases. For Fe I, we 
also used highly-excit ed predicted levels and transitions as rec- 
ommended by Mashonkina et al.l (12010 ). who showed that the 
inclusion of these data is necessary for a realistic representa- 
tion of statistical equilibrium of Fe in the atmospheres of cool 
stars. Photoionization cross-sections for Fe I levels were taken 
from lBau tista ( 19 97|). The detailed description of the Fe model 
is given in Bergemann et al. 

A critical parameter in the statistical equilibrium calculations 
is the efficiency of inelastic collisions with H I. In the absence 
of quantum-mechanical data, we computed the cross-sections 
fo r excitation and ionization by H I atoms from the formulae 
of Id rawi Following the above-mentioned studies, we 

adopted individual scaling factors Sh to the Drawin-type cross- 
sections for each element: Mg (0.05), Al (0.002), Ti (0.05), Cr 
(0), Mn (0.05), Co (0.05). These scaling factors were determined 
by requiring consistent ionization-excitation equilibria of the el- 
ements under restriction of different stellar parameters. For Fe 
and Ba, we used Sh = 0.1, and 0.05, respectively. 

We performed another set of independent non-LTE calcula- 
tions for Na using MAFA GS-OS models and the model atom 
described by iGehren et al.l ([2006), but the mean Na abundance 
only changes by 0.001 dex with respect to the value obtain ed 
using the most up-to-date model atom bv lLind et al.l (1201 lAh . 

Besides the elements above (Li, C, O, Na, Mg, Al, Ti, Cr, 
Mn, Fe, Co, and Ba) for which specific NLTE calculations 
were performed for the present work, we estimate NLTE cor- 
rections for K, Ca, Z n and Zr, usi ng the grid of NL TE cor- 
rections computed b y Takeda et al.l (2002), Mashonki na et al.l 
(l2007h . iTakeda et alJ (12005*) and'Velichko et al. (2010), respec- 
tively. Nevertheless, the differential NLTE abundance correc- 
tions are negligible for these elements (<0.001 dex). 

All the differential NLTE corrections (HIP 56948 - Sun) are 
given in Table [3] except for Li. For this element the differential 
NLTE correction amounts to only -0.001 dex and the adopted 
differential NLTE Li abundance is given in Table |4] 



3.4. The abundance pattern of HIP 56948 

In Fig.|5]we plot the differential abundances [X/H] between HIP 
56948 and the Sun (circles) as a function of equilibrium con- 
densation temperature (Tcond, Lodders, 2003). The fit of [X/H] 
vs. Tcond is shown for volatile (Tcond < 1000 K) and refractory 
(Tcond > 1000 K) elements. □ As can be seen, the element-to- 
element scatter around the fit is extremely small, only 0.004 dex 



' http://www.nist.gov/physlab/data/asd.cfm 

http://kurucz.harvard.edu/ 
" Notice that for the abundance pattern of HIP 56948 the best fit is 
found for a break at T^-onj = 1000 K, while for the average of solar twins 
the break is at T^ond = 1200 K. Although the break point was deter- 
mined by fitting independently the refractory and volatile elements, and 
choosing the break by minimizing the scatter in both sides, with minor 
adjustments to provide the best match at the break point, our results are 
confirmed by a global fitting that assumes two linear functions to fit the 
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H[P 56948 - Sun 
<volatiles> = 0.01 1±0. 002 dex 
<refractories> = 0.020±0.002 dex 
crfvolatiles) = 0.004 dex 
o"(refractories) = 0.007 dex 
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Fig. 5. Abundance pattern of HIP 56948 (circles) versus con- 
densation temperature. The solid line represents the mean abun- 
dance pattern. The error bars are based only on the observa- 
tional uncertainties, which are ~0.005 dex. The low element-to- 
element scatter from the fit for the volatile (cr = 0.004 dex) and 
refractory (cr = 0.007 dex) elements confirms the high precision 
of our work. 
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Fig. 6. Abundance pattern of HIP 56948 (circles) versus con- 
densation temperature. The average of the highly volatile (low 
Tcond) elements C and O is used as reference. The solid line rep- 
resents t he mean ab undance pattern of the 1 1 solar twins studied 
by Mele ndez et al.l (2009) and the dashed line the fit to the abun- 
dance pattern of HIP 56948. Clearly, HIP 56948 is much closer 
to the Sun than to other solar twins. 
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for the volatiles and 0.007 dex for the refractories. Both are of 
the same order as the observational error bars, which are cr„,(obs) 
~ 0.006 dex for volatiles and cr„,(obs) ~ 0.004 dex for refracto- 
ries, hence showing that it is possible to achieve abundances with 
errors as low as ~0.005 dex. If we take into account the errors 
due to uncertainties in the stellar parameters (Table |6), then the 
expected total median errors (including observational errors) are 
0.007 dex for volatiles and 0.008 dex for refractories, i.e., some- 
what higher than the observed scatter (0.004 and 0.007 dex, re- 
spectively) around the mean trends. That means that we may be 
slightly overestimating our errors, and that a more representative 
total error for our abundances is ~0.006 dex. As can be seen in 
Table|6] it is possible to obtain a total error as small as 0.004 dex 
(Si) or 0.005 dex (Ni), and for several elements errors as small 
as 0.006 dex can be achieved (Ca, Cr, Fe, Co, Zn). 

The mean [X/H] ratio of the volatile and refractory elements 
is <[volatiles/H]> = 0.011 (s.e. - 0.002) and < [refractories/H] > 
= 0.020 (s.e. = 0.002), respectively. Thus, the mean difference 
between refractories and volatiles in HIP 56948 amounts to only 
0.009 dex, i.e., HIP 56948 has an abundance pattern very similar 
to solar. 

In Fig.|6]we plot the differential abundances [X/<C,0>] be- 
tween HIP 56948 and the Sun (circles) as a function of con- 
densation temperature. Here the average of carbon and oxygen, 
<C,0>, is chosen as the reference, a s those elements shoul d not 
be depleted in the solar atmosphere ([ M elendez et al.l l2009h . We 
also show with a solid line the mearPi abundance pattern of the 
1 1 solar twins studied by iMelendez et al.l (l2009l) . The dashed 
line represents the mean behavior of HIP 56948 shown in Fig. |5] 
This figure shows clearly that the abundance pattern of HIP 
56948 is much closer to the Sun than to the mean abundance 
pattern of other solar twins. 

The volatile elements with Tcond < 1000 K have a similar 
behavior in HIP 56948 and the Sun, while the refractory ele- 
ments are depleted (with respect to other solar twins) both in HIP 
56948 and the Sun, although they are somewhat less depleted 
(by ~0.01 dex) in HIP 56948. Therefore, it seems that somewhat 
less dust was formed around HIP 56948 than around the Sun. 
Interestingly, in the sample of solar twins studied so far, the Sun 
seems to be the most depleted in refractories. If this peculiar 
abundance pattern is related to the formation of rocky planets 
(Melend ez et al., 2009; Ramirez et al., 2009, 2010; Chambers, 
1201 Ol) . that signature could be used to search for possible candi- 
dates to host terrestrial planets. Below, we discuss in detail vari- 
ous other possibilities that could cause abundance anomalies. 



4. Discussion 

From the determination of isotopic abundances in meteorites it 
has been shown that short-lived radionuclides were present in 
the early stages of the solar system, perhaps due to pollution 
by an asymptotic g iant bra nch (AGB) star (Wasserburg et al., 
Il994t iBusso et al.. 1999; Tri20-Rodn'guez et al., 2009)" or by 
supernova dCameron & Truran. .1977," .Boss & Foster. Il998i 



4. 1 . Abundance anomalies: AGB/SN pollution or terrestrial 
planets? 



lOuellette et aL , 201 Oh . although other causes may be possible. 



Thus, before discussing the terrestrial planet formation scenario, 
we assess if the solar anomalies may be due to other causes such 
as pollution by intermediate or high mass stars. 



whole dataset by assuming that those functions are equal at the break 
point. 

the robust estimator trimean is used 



In 'Mele ndez et al ] ("2009^ we argued against the hypothesis that 
the abundance anomalies found in the Sun could be due to galac- 
tic chemical evolution effects or supernova pollution, as the solar 
chemical peculiarities do not seem to follow those abundance 
patterns. Here, we study in more detail whether the anomalies 
could be due to pollution by an AGB star, thermo-nuclear super- 
novae (SNIa), core-collapse supemovae (SNII) or a hypemova 
(HN), by subtracting the yields of those objects to the solar neb- 
ula (i.e., "de-polluting" the solar abundances) and comparing the 
results to the pattern of the solar twins. 

Following the t e ntative AGB pollution scenario of 
iTrigo-Rodriguez et al.l (l2009l) . we use a dilution factor of 
1 part of AGB material per 300 parts of original solar nebula 
material (equivalent to mixin g 0.0185 Mp of AGB ejecta), for 
which solar abundances from lAnders & Grevesse were 
adopted. Note that the use of the new solar abundances by 
Asplund et al.l (120091) do not have a significant impact on the 
abundance trend, as shown in Fig.|7](top panel). For SNIa, SNII, 
and HN we assume that the same amount of mass as used in the 
AGB scenario is mixed into IM© of solar system material. For 
cons istency with the AGB scena rio we use the solar abundances 
from lAnders & Grevessel (Il989b . 

The amount of material that would actually be injected is 
uncertain. There is no reason that a S N would inject abou t the 
same amount of mass as an AGB star. 'Young et al.l (1201 1 ) note 
that SN are more likely to pollute the entire molecular star form- 
ing cloud, not an individual protosolar nebula. They estimate 
that a ~ 1% enrichment of the protosolar molecular cloud by 
ejecta from SNII can account for the oxygen isotope ratios mea- 
sured in the solar system. However as we are only interested 
in the qualitative impact (trend with condensation temperature) 
that such yields have on a protosolar nebula and not the quanti- 
tative details, the actual amount of pollution from an individual 
SN or AGB model should not strongly impact our conclusions. 
The same can be said for the exact amount of the AGB and SN 
yields, which are inherently quite uncertain. 

A massive AGB star of 6.5 Mq was chosen, as it has a short 
lifet ime (~55 Myr). Our AGB model and yields are d escribed in 
Karakas & Lattanzid (l2007l) . iTrigo-Rodrfguez et al.1 (.2009) and 
iKarakasI The results after removing the AGB pollution 

are shown in Fig. [T] (top panel). As can be seen, contamination 
by a massive AGB star cannot explain the abundance trend with 
Tcond- The main signature from an hypothetical 6.5-M0 AGB star 
would be a large change in nitrogen. 

We have also considered pollution by SN II, HN, and SNIa, 
according to the yields described in Kobavashi et al. (200%^ 
Kobavashi & NomoD (^2009) and .Kobavashi & Nakasatol 
(201 1). Miki et al. (2007) suggest that a supernova of a massive 
star of at least 20 Mq was responsible for the anomalies of 
short-lived radionuclides discussed above, so for the SNII and 
HN, a 25 M0 was adopted. As seen in Fig. |7] the trend with 
condensation temperature cannot be explained by either SNII or 
HN. 

Finally, for the binary system leading to a SNIa, we choose 
a system composed of a SM© of primary star + l.SM© of sec- 
ondary star, with a metallicity of 0.1 Zq for the progenitors, 
which are supposed to have been bom 4.5Gyr ago. The binary 
system evolves to ^IM© C-hO white dwarf plus 1.3 M© of sec- 
ondary star, then to 1.374 Mq SNIa + 0.9 Mq remnant of the 
secondary star Thus, for the contamination by ejecta of SNIa, 
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we include 1.374 of processed metals by SNIa plus (3 + 1.3 
- 0.9 - 1.374) Mq of unprocessed matter in the stellar winds. 

Again, pollution by SNIa cannot explain the peculiar solar 
abundance pattern (Fig.|7ll. We also tried other combinations for 
the progenitors of SNIa, but these yields also do not reproduce 
the trend with condensation temperature shown in Fig. |2] 

Therefore, we conclude that the peculiar solar abundance 
pattern cannot be due to contamination by AGB stars, SNIa, 
SNII or HN, as predicted by state-of-the-art nucleosynthesis 
models. We emphasize that our results do not rule out that a 
SN or AGB star contaminated the proto-solar nebula to cause 
the observed isotopic anomalies in meteorites, rather, we dis- 
card pollution as a viable explanation of the peculiar elemental 
abundances in the Sun. 



4.2. The abundance pattern of HIP 56948 and terrestrial 
planets 

As discussed above, pollution from stellar ejecta cannot ex- 
plain the anomalous solar abundance pattern. A possible ex- 
planation for the peculiarit ies is the formation of terres- 
trial planets jM elendez et a ll. l2009HRanurez et all l2009l 1201 Ot 
iGustafsson et al.1. l2010t iChambersll2010h . The same planet for- 
mation scenario may be applied to HIP 56948. To verify this, 
we obtained th e abundance ratios of the 11 solar twins of 
iMelendez et alJ (I2b09i) relative to HIP 56948, again using the 
average of C and O as reference. The resulting mean trend is 
represented by a dashed line in Fig. |8] The trend of solar twins 
relative to the Sun is shown by a solid line. 

Thus, the same teiTestrial planet formation scenario could 
be applied for HIP 56948, except that the Sun has -0.01 dex 
smaller refractory-to-volatile ratio than HIP 56948. So, the Sun 
could have formed more rocky planets than HIP 56948, or per- 
haps slightly more massive rocky planets than in HIP 56948. In 
any case, the overall amount of rocky material may have been 
higher around the Sun than around HIP 56948. 

Recently, Chambers (2010) has shown that a mix of 4 Earth 
masses of Earth-like and meteoritic-like material, provides an 
excellent element-to-element fit for the solar abundance anoma- 
lies for about two dozen chemical elements^ Interestingly, 
the above mixture reproduces the anomalies significantly bet- 
ter than a composition based either only on Earth material or 
only on carbonaceous-chondrite meteorites. Using his detailed 
abundance pattern of Eaith and carbonaceous chondrites, we can 
check whether the same mixture could fit HIP 56948. Indeed, 
as shown in Fig. [8] a 3-Earth-mass mixture of Earth plus mete- 
orites (filled circles) provides an excellent fit to the abundance 
pattern of the solar twins relative to HIP 56948 (dashed line). 
The element-to-element scatter of the Earth/chondrite mix with 
respect to the fit (dashed line) is only 0.005 dex. Notice that 
the chemical elements included in Fig. [8] are not necessarily the 
same ones studied in HIP 56948 because for the Earth/meteoritic 
mix we can includ e other elements such as fluorine and uranium 
(IChambersll20T(ih . 

The abundance pattern of the solar twins with respect to the 



Sun (Fig. [8] solid line ), seems to requir e more than the four Earth 
masses suggested by IChambersI (I201Q). No te that we are nor- 
malizing to <C,0>, while 



Chambersl(l20"Tol) normalizes to iron. 



Notice that although this indicates that a certain amount of both 
Earth-like and meteoritic-like material may have been removed from 
the Sun (in comparison with the solar twins), this does not imply nec- 
essarily that the removed material was employed to form the terrestrial 
planets and asteroids. 



and th at the mean solar abundance trend adopted by IChamber"sl 
(20 10*) may be slightly different than ours, hence some differ- 
ences are expected. The fit for the Sun (cr = 0.010 dex) is not as 
good as for HIP 56948 (cr = 0.005 dex). In total 6 of rocky 
material is needed to fit the Sun, while 3 Me is needed to fit HIP 
56948. So, perhaps about twice more rocky material was formed 
around the Sun than around HIP 56948. 

Interestingly, both the Sun and HIP 56948 seem to require 
about the same mass of Earth-like material (2 and 1 .5 M®, re- 
spectively), but the amount of CM chondrite material seems 
much higher for the Sun (4 M^) than for HIP 56948 (1.5 Me). 
Although 4 Earth masses of meteoritic material may seem too 
high considering the current mass of the asteroid belt, var- 
ious models predict that the belt has lost most of its mass 
(e.g.,'Weidenschining', '1977'; 'Wetherill', '1989", '1992; 'Petit et aU 
1999; Chambers & Wetherill, 2001; Minton & Malhotra. 201^ 
In particular, the simulations by Chambers & Wetherilll (120011) 
show that a few Earth masses of material from the belt can be 
removed, either via collisions with the Sun or ejected from the 
solar system. In contrast, material in the terrestrial-planet region 
had a good chance of surviving (Chambers & Wetherill, 200l!). 
Thus, albeit a large mass of material may have formed in the 
asteroid belt, most of the initial mass has already been removed. 

The 2 Me of Earth-like material seen in the chemical com- 
position of the Sun is comparable with the total mass of the ter- 
restrial planets in the solar system, and only somewhat higher 
than the amount of Earth-like material around HIP 56948. As 
shown above, the different slopes of the abundances of volatile 
and refractory elements (with condensation temperature) may 
be used to constraint the type of solid bodies that were origi- 
nated as a result of planet formation using the different abun - 
dance patterns of the Earth and meteorites (IChambersi l2010l) . 
Thus, the careful analysis of stellar chemical compositions offer 
the thrilling prospect of determining which type of rocky objects 
were formed around stars. 

In line with our above findings on the similarities between 
HIP56948 and the Sun, the radial velocity monitoring of HIP 
56948 (Section 2.2) shows no indication of inner (< 3 A.U.) gi- 
ant planets, as massive as Saturn or Jupiter. Thus, the inner re- 
gion around HIP 56948 can potentially host terrestrial planets. 
The remarkable chemical similarities between HIP 56948 and 
the Sun, also suggest that HIP 56948 may be capable of hosting 
rocky planets. 

Also, metal-rich solar analogs without close-in giant planets 
seem to have an abunda nce pattern closer to solar than stars with 
detected giant planets ( Mel6ndez et al.Ll200 9)Fl 

Our findings open the truly fascinating possibility of identi- 
fying Earth-mass planets around other stars based on a careful 
high resolution spectroscopic analysis of stellar chemical com- 
positions. Once the Kepler mission (e.g. iBorucki et al.L l2010h 
announces the discovery of Earth-sized planets in the habitable 
zones of G-type dwarfs, our planet signature could be verified 
by high precision chemical abundance analyses of those stars. 
Although they are relatively faint (confirmed Kepler planet- 
hosting stars have a mean magnitude of /Tp - 13.8+1.5), recently 
On ehag et al.l (1201 ll) have shown that it is possible to achieve 
high precision (0.03 dex) differential abundances even in stars 
with V=15. They analyzed the faint (V = 14.6) solar twin M67- 
1 194 using VLT/FLAMES-UVES, showing that it has a remark- 
able chemical similarity to the Sun. 



N otice that, as discussed in lSchuler et al.l ll201lh and lRamfrez et al.l 
iZOlOh . the interpretation of abundance trends for metal-rich stars is 
complicated by Galactic chemical evolution processes. 
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4.3. The mass, age, luminosity and radius of HIP 56948 

Contrary to commonly thought, reasonable estimates of the ages 
of main-sequence stars can be obtained using standard isochrone 
fitting techniques{3 provided the isochrones are accurate (i.e., 
calibrated to reproduce the solar age and mass) and the stellar 
parameters Teff, \ogg, [Fe/H] are known with extreme precision, 

as illustrated in Fig. |9] 

W e used a fine grid of Yonsei-Yale isochrones ( Yi et al.L 
I2OO I: Kimetal., 2002; Demarque et al., 2004) with a step 
A[Fe/H] = 0.01 dex around solar metalHcity (-0.15 < [Fe/H] 

< +0.15) and a step of A[Fe/H] = 0.02 dex elsewhere. We adopt 
[a/Fe] = for [Fe/H] > 0, [a/Fe] = -0.3 x[Fe/H] for -1 < [Fe/H] 

< 0, and [a/Fe] = +0.3 for [Fe/H] < -1. The isochrones include a 
dependence between helium abundance Y an d metallicity Z with 
a slope of 2, Y = 0.23 + 2 Z (lYi et al.Ll200l . 

Our grid has been normalized to reproduce the solar age and 
mass from the input solar parameters (Teff = 5777 K, log g - 
4.44, [Fe/H] - 0.0). The normalization factor was found by per- 
forming small offsets around the solar Tgif, log g and [Fe/H] to 
see which offsets better reproduce the solar age and mass. We 
found that shifts in Teff and log g are not needed, as the best 
compromise solution can be found with only a small shift of 
-0.04 dex in the observed [Fe/H], meaning that the models are 
off by +0.04 dex in metallicity. Thus, the input metallicity used 
to compare with the isochrones is: 



[Fe/H] 



input 



[Fe/H] - 0.04 dex. 



(5) 



That normalization gives a mean solar mass of 1.000 (to 
within 0.003 Mq) and a mean solar age of ~4.5 Gyr (to within 
0.2 Gyr), which is in exce llent agreemen t with t he age of the 
solar system (-4.567 Gyr, IConnellv et al.[ l2008t lAmelin et al.L 
I2OIQ) . The calibration of the models to the solar mass and age 
is valid for a broad range of errors of 10-140 K in T^g and 0.01- 
0.10 dex both in log g and [Fe/H] . Thus, after the zero-point shift 
of Eq.|5]is applied, our resulting masses and ages are accurate. If 
the input errors are much higher than those indicated above, the 
zero point of the solutions need to be revised^ 

The isochrone points are characterized by effective temper- 
ature (T), logarithm of surface gravity (G), and metallicity (M), 
with a step in metallicity of 0.01 dex around [Fe/H] - 0. 

We obtained an estimate of the age of HIP 56948 from its 
isochrone age probability distribution (APD): 



dPiags) = — ^ y p (Feff, logg, [Fe/H], T, G, M) 
A(age) ^ 



(6) 



A(age) 



where Teff, log [Fe/H] are the observed stellar parameters, 
A(age) is an adopted step in age from the grid of isochrones, 
and: 



p cc exp[-(reff - Tf/2iAT,fff] x 
cxp[-(logg-Gf/2(Aloggf]x 
exp[-([Fe/H] - Mf /2iA[Fs/U]f] 



(7) 



Baumann et"^ 12 010): Melendez et al. I2010Q: Bens bv et alJ 
RamfrezeTaLriTOl 1.) : .Chaname & Ramirez. t201Z) . for difFer- 



se e 
( l20TTh : 

ent applications of our isochrone ages. 

For errors of 175 K in Teff and 0.1 dex both in log g and [Fe/H], 
the solar mass and age would be about 0.99 Mq and 4.8 Gyr, while for 
an error of 250 K in T^ff and 0.15 dex both in log g and [Fe/H], the 
solar mass and age would be about 0.98 Mq and 5.5 Gyr, respectively. 
However, for typical errors of abundance analysis of cr(T^ff ) < 150 K, 
cr(log g) < 0. 1 dex and (T([Fe/H]) < 0.1 dex, no zero-point corrections 
(besides that already used in Eq.|5} would be needed to obtain accurate 
masses and ages. 



The errors in observed stellar parameters are AT^fi, etc. The sum 
in Eq.|6]is made over a range of isochrone ages and in principle 
all values of T,G, M. In practice, however, the contribution to the 
sum from isochrone points farther away than T^ff + 3Aretf, etc., 
is negligible. Therefore, the sum is limited to isochrone points 
within a radius of three times the errors around the observed 
stellar parameters. A similar formalism allows us to infer the 
stellar mass. The probability distributions are normalized so that 
YjdP - I. The most probable age and mass are obtained from 
the peaks of these distributions while 1 cr and 2 cr Gaussian-like 
lower and upper limits can be derived from the shape of the prob- 
ability distributions. 

Fig. [H shows the APDs of the Sun and HIP 56948. For HIP 
56948 we adopted errors from our differential analysis relative to 
the Sun of 7 K in Teff, 0.02 dex in \ogg, and 0.01 dex in [Fe/H], 
while for the Sun we adopted errors of 5 K in T^ff and 0.005 dex 
in \ogg and [Fe/H]. Although the errors for the Sun are overes- 
timated, they allow us to obtain a smooth APD. 

Using our precisely determined stellar parameters for 
HIP 56948, we derive a most probable age of 3.45 Gyr The 1 cr 
range of ages is 2.26^.12 Gyr whereas the 2 cr range of ages is 
1.25^. 92 Gyr Notice that although the adopted He abundance 
may have some influence on the derived age, our small error bars 
in the stellar parameters of HIP 56948, rule out radically differ- 
ent He abundances. For relatively small changes in He, the effect 
on the derived stellar age is relatively minor This is discussed in 
detail in appendix D, where stellar tracks with different He abun- 
dances are presented. The mass is much better constrained; for 
HIP 56948 we derive 1 .020 + 0.016 Mq (2 a error). 

In Fig. [To] we also show the APD of HIP 56948 assuming 
that its parameters were obtained using standard methods, for 
example obtaining T^ff from photometric data and log g from the 
Hipparcos parallax. The standard method imply an error of at 
least 100 K in T^^ due to uncertainties in the zero-point ofcolor- 
Teff relations (Casagra nde et al.L 1201 Ot iMelendez etall. 1201 OBh . 
while the error in Hipparcos parallax and typical errors in mass, 
Teff and V magnitude imply a total error of about 0.07 dex in 
log g. For the error in [Fe/H] of a typical abundance analysis 
we adopt 0.05 dex. Clearly, in this case the age of HIP 56948 is 
not well constrained. At most, we can say that the star is likely 
younger than about 8 Gyr Thus, another advantage of studying 
solar twins using very high quality spectroscopic data and strict 
differential analysis is that useful estimates of their ages can be 
obtained with the classical isochrone method, which can be a 
valuable asset for a variety of studies (e.g., Baumann et al. 2010). 

The luminosity of HIP 56948 can be estimated from: 

log L/Lo = log(M/MQ) - (logg-4.44) + 4 log(reff/5777)(8) 

where M is the stellar mass and L the luminosity. Using our pre- 
cise stellar parameters and mass we find L - 0.986 ± 0.05 1 Lq. 
Thus, the luminosity of HIP 56948 is essentially solar and there- 
fore the extent of its habitable zo ne should be similar to the 
habitable region around the Sun (Kastin g et al.L Il993h . From 
L = AnR^o-Tl^, we find a radius of ^ = 0.987 ± 0.023 i^o, i.e., 
solar within the error bars. 



4.4. Further constraints on the age of HIP 56948 

Additional insight on the age of HIP 56948 can be obtained from 
its chro mospheric activity (Soderblom, 2010 ), lithium abun- 
dance (do Nascimento et al., 2009t iBaumann et all 120101) and 
gyrochronology (Barnes, 2007). 

Determination of stellar ages based on chromospheric activ- 
ity may only be valid for solar type stars younger than ~2 Gyr 
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(iPace & PasquiniL l2004t iPace et all 120091: IZhao et all l201lh . 
as older stars show a low activity level which changes lit- 
tle with increasing age. Chromospheric activity is thus mainly 
useful to distinguish young and old stars. A measurement 
of chromospheric ac ti vity f or HIP 56948 was obtained by 
iMelendez & RamfrezI (l2007h based on observations taken in 
April 2007, resulting in a low chromospheric 5 -value of S - 
0.165 (±0.013) in the Mount Wilson scale. Using another ob- 
servation taken in November 2007 we find a similar value, S 
= 0.170 (+0.013). Both values are as low as the the mean 
chromospheric activi ty index in the Sun, < Sq > - 0.179 
dBaliunas et aUll995h . The low 5 -value of HIP 56948 suggests 
that it should be older than ~2 Gyr. 

Although it is well-kn own that lithium st eeply decays with 
age in very young stars (ISoderbloml 1201 Ob . only recently it 
has been o bservationally shown th a t the decay continue s for 
older ages (IMelendez et al.L l2010At iBaumann et"an [2010), as 
already predicted by several models of non-stan dard Li depletion 
(iMontalban & Schatzman. 2000; Charbonnel & TalonL 12005 



Age 



IS 



a key parameter for SETl programs 



do Nasci mento et al.i , i2009i : iXiong & Dena , i2009i : iDenissenkov 



2010). Unfortunately we cannot determine stellar ages with the 
required precision to discern to what extent there is a spread 
(or not) in the age-lithium re lation, although the most pre- 
cise values available (this work: iRamfrez et all |20 1 ll Melendez 
et al. 2012, in preparation) show relatively little dispersion 
(Fig. fTTI ). Theoretical Li-age relations can be used to check 
if a lithium age is compatible with the age determined from 
isochrones. In Fig. [TT] we compare s everal Li theoretical 
tracks jCh arbonnel & Talon, 2005; do Na scimento et al.L 120091: 
IXiong & Dengt l2009l: iDenissenkovt 12010" ) with our NLTE Li 
abundances and isochrone ages. As can be seen, the agreement 
is excellent. If we were to derive an age for HIP 56948 based on 
the theoretical Li-age relations, it would be 3.62±0.19 Gyr, i.e., 
almost the same age determined using isochrones. This compar- 
ison give us further confidence that HIP 56948 is about 1 Gyr 
younger than the Sun. Although the scatter of the Li age ob- 
tained using different Li tracks is relatively small (0.19 Gyr), we 
conservatively assign an error of 1 .0 Gyr to the Li-age, to take 
into account any possible observational spread of the Li-age re- 
lation around solar age (Fig. fTTI. 

Stellar rotation can be used to estimate a rotational age 
dBarnesl l2007t ISoderblomi |2010|) . Unfortunately there is no in- 
formation on the rotation period of HIP 56948 yet. Nevertheless, 

V sin / can give us an upper limit on the rotation period, thus al- 
lowing to infer an upper limit on the age. The determination of 

V sin / is based on the differential line broadening between HIP 
56948 and the Sun, as described in the appendix E. 

We infer for HIP 56948 v sin //vsin /q = 1.006 ± 0.014, or A 

V sin i = -hO.013 + 0.026 km s"' (or +0.032 km s"' including the 
error in macroturbulence), i.e., HIP 56948 seems to have about 
the same rotation velocity as the Sun. Within the uncertainties 
we infer that HIP 56948 cannot be older than th e Sun. Using the 
relation between rotation period and age given in'Me lendez et alj 
(12006) and .Barnes, (.2007,) . we find an upper limit of age < 4.7 
Gyr. 

Fig. [121 summarizes our findings on the age of HIP 56948. 
Our precise stellar parameters and differential isochrone analy- 
sis result in an age of 3.45 ± 0.93 Gyr (1-cr error). The somewhat 
higher Li abundance of HIP 56948 with respect to the Sun in- 
dicates an age of ~3.62+1.00 Gyr. Chromospheric activity gives 
a lower limit of 2 Gyr, while v sin / suggests an upper limit of 
4.7 Gyr. In Fig. [12] we show the combined age probability dis- 
tribution. Based on all the above indicators we suggest an age of 
3.52+0.68 Gyr for HIP 56948. 



(iTurnbuU & Tarten l2003l) . as stars only 1-2 Gyr old may 
not have had enough time to develop complex life. Life on Earth 
apparently appeared within the first billion year of the Earth's 
formation ('Schopf, '1993"; 'Moizsi s et all 119961: iMcKeegan et all 
2007; Abramov & Moizsis, 2009), but there is no consensus on 
the exact date. Undisputed evidence for life can be traced back 
to about 2.7 Gyrs ago (e.g., see review by Lopez-Ga rcia et all 
200 6|). Yet, c omple x life only appe a red ab o ut 0.5-1 Gyr ago 
(Wrav et al.', '1996^ 'S eilacher etall Il998t iRasmussen elall 
2002; Marshall, 2006). HIP 56948 is about one bilUon year 
younger than the Sun, so assuming a similar evolution path as 
that of life on Eaith, complex life may be just developing (or 
already sprung if complexity elsewhere can arise earlier than on 
Earth) in any hypothetical Earth that this star may host 

5. Conclusions 

We have shown that using spectra of superb quality coupled to a 
fully differential analysis of solar twins and the Sun, it is possible 
to achieve measurements errors as low as 0.003 dex for several 
elements (Si, Ca, Ti, Cr, Fe, Co, Ni). Considering also the uncer- 
tainties in stellar parameters, we achieve an unprecedented accu- 
racy of only ~0.005 dex (1 %) in relative abundances for some 
elements and < 0.01 dex for most elements. This is almost one 
order of magnitude better than state-of-the-art works in terms o f 
absolute abundances (e.g.. lAsplundll20()5tlAsplund et am2009l) . 

The star HIP 56948 is remarkably similar to the Sun in many 
different aspects. The effective temperature, log g, metallicity 
and micro turbulence are very similar. The similarities also ex- 
tend to its detailed chemical abundance pattern. The volatile el- 
ements are in excellent agreement, but the refractory elements 
are slightly (0.01 dex) more enhanced relative to the volatile 
elements in HIP 56948. From the comparison with the abun- 
dances of Earth-like and chondrite-like material, we infer that 
about twice as much rocky material may have formed around 
the Sun than around HIP 56948, albeit the amount of Earth-like 
material is comparable for both stars (~2 Me). The mass, lu- 
minosity and radius of HIP 56948 are essentially solar within 
the uncertainties. Lithium is severely depleted in HIP 56948, but 
not as much as in the Sun, as expected for a solar twin some- 
what younger than the Sun. Finally, our precise radial velocity 
data shows that the inner region around HIP 56948 is free from 
giant planets, making thus more likely the existence of teiTes- 
trial planets around this remarkable solar twin. Considering its 
similarities to our Sun and its mature age, we urge the commu- 
nity to closely monitor HIP 56948 for planet and SETI searches, 
and to use other techniques that could further our knowledge 
about HIP 56948, such as for example asteroseismology, that 
have provided imp ortant constraints for the solar twin 18 Sco 
(iBazotet all [201 ll) . 

The abundances anomalies we have discussed here cannot 
be explained by contamination from AGB stars, SNIa, SNII or 
HN (Sect. 4), or by Galactic chemical evolution process es or 
age effects (IMelendez et all l2009l) . iKiselman et alJ ilQl ih have 
shown that the peculiar abundance pattern cannot be attributed 
to line-of-sight inclination effects. Also, the abundance trend do 
not arise due to the particular reflection properties of asteroids 
(appendix B). Although the abundance peculiarities may indi- 
cate that the Sun w as born in a massive open cluster like M67 
(lOnehag et al.L l201 Ih . this explanation is based on the analysis 
of only one solar twin. The Uppsala group is leading a high pre- 
cision abundance study of other solar twins in M67, in order to 
confirm or reject this hypothesis. So far the best explanation for 
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the abundance trend seems to be the formation of terrestrial plan- 
ets. The Kepler mission should detect the first Earth-sized plan- 
ets in the habitable zones of solar type stars. We look forward to 
use 8-lOm telescopes to perform careful differential abundance 
analyses of those stars, in order to verify if our chemical signa- 
tures indeed imply rocky planets. 
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Fig. 7. Abundance ratios obtained after de-polluting the so- 
lar nebula from contamination by an AGB star (circles), 
SNll (stars), hypernova (squares) and SNla (triangles). In the 
top panel it is shown the effect of adopting differe nt so- 
lar abu ndances (open c i rcles: [Anders & Grevessel (Il989b : filled 
circles: 'Asplund et al.J (l2009l) V The solid line represents the 
mean abundance pattern of 1 1 solar twins relative to the Sun 
jMelendez et al.L 12009). None of the pollution scenarios can ex- 
plain the trend with condensation temperature. For clarity, the 
abundance ratios of SNlI, HN and SNla have been divided by 4, 
3.5 and 14, respectively. The chemical elements that change the 
most are labeled. 
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Fig. 8. Composition of the solar twins with respect to HIP 56948 
(dashed line) and the Sun (solid line). The open circles show 
the effect of a dding 6 Mm of a m ix of Earth-like and meteoritic- 
like material ( Chambers^ 1201 Ol) to the convection zone of the 
present Sun dAsplund et aL, 2009), and the filled circles the ef- 
fect of adding 3 of rocky material to the convection zone of 
HIP 56948. Abundances are normalized with respect to <C,0>. 
Both the Sun and HIP 56948 require ~2 Me of Eaith-like mate- 
rial, while the Sun requires much larger quantities of chondrite 
material. Notice the very small (<0.01 dex) element-to-element 
scatter. 
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Fig. 9. Location of the Sun and HIP 56948 on the HR dia- 
gram. Note the very small range of stellar parameters. Solar- 
metallicity isochrones of 1.5, 3.3, 4.6, and 5.5 Gyr are shown 
(dotted lines). The high precision of our derived stellar parame- 
ters for HIP 56948 allows us to infer a reasonable estimate of 
its age from the theoretical isochrones, even though they are 
densely packed in this main-sequence region. 
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Fig. 10. Age probability distributions (APDs) for the Sun (dotted 
line) and HIP 56948 (solid line). The dashed line corresponds to 
the APD that we would obtain for HIP 56948 if we had derived 
its stellar parameters from photometric data and Hipparcos par- 
allax (the "standard" method) instead of performing our very 
precise strict differential analysis. Clearly, our differential ap- 
proach gives much better results. 
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Fig. 11. Li vs. age for the Sun (0) and HIP 56948 (square) based 
on our NLTE Li abundances and isochrone ages, and for 16 Cyg 
B (filled circle) and 18 Sco (triangle, 2 .7 Gyr^^]^^^), based on simi- 
lar quantities bv lRamfrez et al.l (1201 ih and Melendez et al. (20 12, 
in preparation). The total error bar (+cr) of the Li abundance is 
about the size of the symbols, while the error bars in age are 
show n by horizontal lines. For comparison we show the mod - 
els by'Charbonnel & Talon' (2005); 'do NascjrnentoeLalJ(l2009t); 
IXiong & Deng (2009); Denissenkov (2010), shifted in Li abun- 
dance by 0.00, -0.03, -0.15, -0.05 dex, respectively, to repro- 
duce our observed NLTE solar Li abundance. The age of HIP 
56948 based on Li tracks is in perfect agreement with the age 
obtained from isochrones. 



15 



Jorge Melendez et al.: The best solar twin 



1.0 



0.8 



0.2 



0.0 



' III:-^ 56948 








isochrones 








1 i Lh i u rn 




ago = 3.52 + 0.68 Gyr 




combined 


1 




- 
- 






S \ 




- 


activity 

- 


^/ / 
7 

■f 


■a\ 

'■.\\ 
V 






..•>■ 

■■/ 


•V 

: / 

/ 


V sin 


\ 

. \ 
■. \ 








■ ■.. \ 















012345678 



Age (Gyr) 

Fig. 12. Age probability distributions for HIP 56948 based on 
isochrones (dotted line) and lithium abundance (dashed line). 
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age probability distribution (solid line) is centered at 3.52 Gyr 
and has cr=0.68 Gyr. 
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Appendix A: Is the Tcond-abundance trend real ? 

We are just starting the era of high precision (0.01 dex) abun- 
dances studies, therefore the casual reader may question how 
real or universal is the abundance trend (<solar twins> - Sun) 
with condensat i on tem perature. The trend was first found by 
iMelendez et al] ( l2009h . who determined a Spearman correlation 
coefficient of rg = +0. 91 and a negligible probability of only ~ 
10"' of this trend to happen by pure chance. These results, based 
on Southern solar twins observed at the Magellan telescope, are 
reproduced in the left-upper panel of Fig. lA.ll where the average 
abundance ratios of the solar twins is plotted against condensa- 
tion temperature. Additional independent works are also shown 
in Figure lA.ll where a line representing the mean trend found 
by Melend ez et al. (2009), is superimposed upon the different 
samples. 

The independent study by iRamfrez et al.l (l2009h . using 
McDonald data of Northern solar twins, follows the same trend 
(Fig. I A. II upper-right pan el), as well as the average of six 



2003 



Allende Prieto et al 



indep endent s amples ( Reddv et al. 

l2004t iTakedaL 120071: iNeves et all l2009t iGonzalez et all 120 IC 
Bensby et al., i n prep aration) of solar analogs in the literature 
(Ramirez et al.Ll20Tot) . as shown in the left-middle panel of Fig. 
ED 

The revision and extension (iGonzalez Hernandez et al.L 



1201 Ol) of the abundance analysis by Neves et al. ( 20091) of the 
HARPS high precision planet survey, also follows the same 
trend, except for a minor global shift of only -0.004 dex, as 
illustrated in the right-middle panel of Fig. lA.ll In this panel we 
show the average abundance ratio^ of 15 HARPS solar twins 
with Jeff, log g and [Fe/H] within ±100 K, ±0.1 dex and +0.1 
dex of the Sun's stellar parameters. The agreement between the 
solar twin patt ern of flonzalez Hernandez et al. (2010) and the 
mean trend of IMelendez et al.l (l2009l) is good, except for the 
element O and to a lesser extent for S, for which it is diffi- 
cult to determine precise abundan ces. In particular, notice that 
IGonzalez Hernandez et al. I (l20T0l) derived oxygen abundances 
from the [OI] 630nm line, which is badly blended with Nil. Also, 
notice that the [OI] feature is the weakest employed by them. 
Since their oxygen abundance is based on a single line, which 
is the weakest of all features analyzed by them, and that [OI] is 
blen ded with Nil, it is n atural to expect that the O abundances in 
IGonzalez Hernandez et al. (2010) have the largest uncertainties. 
Furthermore, the HARPS spectra taken for planet hunting shows 
some contamination from the calibration arc, so that abundances 
derived from only one single feature should be taken with care. 

Regarding the analysis of individual stars, most of them 
also display the abundance trend, as show n for example for 
four solar twins in Fig. 2 o f Ramirez et (l2009b and 1 1 so- 
lar twins in Fig. 4 of Melendez e t al.i (120091 ). Two new examples 
are shown in Fig. lA.ll In the bottom-left panel we shown the 
average abundance of t he pair of solar analogs 16 Cyg A and 
B (.Ramirez et al.Ll201 1 ). based on high resolution (R = 60,000) 
and high S/N (~400) McDonald observations. As can be seen, 
this pair also follows the abundance trend, after a minor shift 
of -0.015 dex. In the bottom-right panel we show the abun- 
dance ratios of the solar twin 18 Sco (Melendez et al. 2012, in 
preparation), based on high quality (R = 110,000, S/N ~ 800) 
UVES/VLT data. It is clear that the abundance trend is also fol- 
lowed by 18 Sco, after a shift of only -1-0.014 dex in the abun- 
dance ratios. Similar results are obtained using HIRES/Keck 
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" The individual abundance ratios with errors larger than .1 dex 
were discarded from the data of lGonzalez Hernandez et all (1 201 Oh when 
computing the average [X/Fe] values. 



Fig.A.l. [X/Fe] ratios (from carbon to zinc) vs. condensa- 
tion temperature for different samples of solar twins and so- 
lar analogs. The solid line represents the mean trend found by 
Mele ndez et al., (2009). upper-left (fil l ed circ les): Southern sam- 
ple of solar twins by IMelendez et all (l2009l): upper-right (ope n 
circles): Northern solar twin sample by Ramirez et al.l (1200 9l) : 
middle-left (squares): average of six diffe rent literature sam- 
ples of solar analogs ( Ramfrez et al.L |20I0|) : middle-right (pen- 
tagons): average of 15 so l ar twi ns in the HARPS sample of 
IGonzalez Hernandez et all (1201 Ol) . after a shift of -0.004 dex; 
lower-left (triangles): average of the pair of solar analogs 16 
Cyg A and B (Ramfrez et al., 201 1), after a shift of -0.015 dex; 
lower- right (stars): abundance pattern of the solar twin 18 Sco 
(Melendez et al. 2012, in preparation), after a shift of -1-0.014 
dex;. 



data (Melendez et al. 2012, in preparation; see also appendix 
B). 

Thus, all recent high precision abundance studies based on 
different samples of solar twins and solar analogs in the Southern 
and Northern skies, using different instrumentation (Tull Coude 
Spectrograph at McDonald, MIKE at Magellan, HARPS at La 
Silla, UVES at the VLT, HIRES at Keck), all show the abun- 
dance trend. In concl usion, it seems that the reali ty of the abun- 
dance trend found by 'Mele ndez et al.l (l2009l) and iRarmrez et all 
(.2009,) . is well estabhshed. 



Appendix B: Test of our precision using the 
asteroids Juno and Ceres 

The referee suggested that we test our method using observa- 
tions of two asteroids of different properties, obtained with the 
same instrument and setup, in order to show whether our very 
small standard errors (~0.005 dex) are adequate to estimate the 
observational uncertainties, as well as to look for potential sys- 
tematic problems with the asteroid Ceres. Although we have not 
acquired such data yet, we do have observations of the solar twin 
18 Sco and two different asteroids observed with different in- 
struments: high quality UVES spectra of 18 Sco and the asteroid 
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Fig. B.l. Differences between the solar abundances obtained 
with the asteroids Juno and Ceres, obtained through (Juno - 
Ceres) = [X/H](18 Sco - Ceres) - [X/H](18 Sco - Juno). The 
sohd line shows the mean difference and the dotted lines show 
the element-to-element scatter. 

Juno (R = 1 10,000 and S/N ~ 800) and high quality HIRES spec- 
tra of 18 Sco and the asteroid Ceres (R = 100,000, S/N ~ 400). 
Juno is a S-type a steroid and Ceres is a C-type asteroid (e.g., 
iDeMeo et al.L l2009i) . therefore they have very different spectral 
properties and the relative analysis of 18 Sco to both Juno and 
Ceres should reveal if there is any problem in using their re- 
flected solar light. 

The analysis has been performed as described in Sect. 3. Our 
preliminary LTE results for the stellar parameters of 18 Sco are 
Teif = 5831 + 10 K, log g = 4.46+0.02 dex, [Fe/H] = 0.06+0.01 
dex. Full details of the abundance analysis will be published 
elsewhere. In Fig lB.ll we show the difference between the [X/H] 
ratios obtained in 18 Sco using the Ceres and Juno asteroids, 
[X/H]i8Sco-Ceies " [X/H] issco-juno, or in Other words the abun- 
dance difference (Juno - Ceres). Notice that the same set of lines 
was used for both analyses. The error bars shown in Fig lB.l l are 
the combined error bar based on the standard error (s.e.) of each 

analysis, i.e, error = ^ls.e.\^^^^_^^^^ + s.e.^gsco-juno- 

As can be seen, the standard errors fully explain the small 
deviations of the (Juno - Ceres) abundance ratios. The mean dif- 
ference <Juno - Ceres> is only 0.0017 dex, and the element- 
to-element scatter is only 0.0052 dex, meaning that each of the 
individual analyses should have typical errors of about 0.003 - 
0.004 dex. The agreement is very satisfactory considering the 
different instrumentation employed and that the comparison be- 
tween Juno and Ceres is done through a third object (the solar 
twin 18 Sco). Also, notice that there is no meaningful trend with 
condensation temperature. The test performed here strongly sup- 
ports for our high precision and removes the possibility that the 
abundance trend may arise due to the particular properties of as- 
teroids. 

Besides the potential applications of high precision differ- 
ential abundance techniques to study the star-planet connec- 



tion, _diesejechniguesare^^ giving new insights in other ar- 
eas. Nissen & Schusteil (120101) achieved uncertainties of 0.03 
dex in [Mg/Fe] and only 0.02 dex in both [Ca/Fe] and [Ti/Fe], 
showing a clear separation of the halo into two distinct popu- 
lations with different [g/Fe ] ratios. Regarding globular clusters, 
iMelendez & CohenI (l2009h have shown that CN-weak giants in 
M71 show a star-to-star scatter in [O/Fe] and [Ni/Fe] of only 
0.018 dex, while [Mg/Fe] and [La/Fe] have a scatter of 0.015 
dex. The star-to-star scatter in metallicity ([Fe/H]) is only 0.025 
dex. Even a lower star-to-star scatter is found among "globu- 
lar cluster star twins" (stars within +100 K of a globular cluster 
standard star) of NGC 6752. Using superb spectra (R - 1 10,000; 
S/N = 500) obtained with UVES on the VLT ( Yong et al., 2001 
l2005h and applying similar techniques to those presented in Sect. 
3, Yong et al. (2012, in preparation) have found an unprecedent- 
edly low star-to-star scatter of only 0.003 dex in the iron abun- 
dances among NGC 6752 star twins, reveaUng chemical homo- 
geneity in this cluster at the 0.7% level. 



Appendix C: Determination of stellar parameters 

As mentioned in Sect. 3, the excitation and ionization equilib- 
rium do not depend only on Tetf and log g, respectively. There 
is some dependence with other stellar parameters, but to a much 
lesser extent, such that a "unique" solution can easily be obtained 
after a few iterations. In practice, considering the weak degen- 
eracies, a first guess of the effective temperature can be obtained 
computing the slope at three different Teff (e.g., in steps of 50 K 
and at fixed solar log g) at the best microturbulence velocity (at 
a given Teff and log g). Then a linear fit is performed to Teff vs. 
slope (see Fig. IC.lb to find the effective temperature at slope = 
0. Then, for this T^-s we can run three models with different log 
g (e.g., in steps of 0.05 dex in log g) in order to find the best 
surface gravity, by fitting log g vs. A''"' (see Fig. IC.2I ). and for 
each model the microturbulence is obtained. This leads to the 
first guess of Teff, log g and v,. Further iterations at smaller steps 
(down to 1 K in T^s , 0.01 dex in log g and 0.01 km s"' in v,) can 
quickly lead to the best solution that simultaneously satisfy the 
conditions of differential spectroscopic equilibrium (eqs. 2-4). 

In Fig. IC.ll we show that the excitation equilibrium pro- 
vides a precise Teff . In this figure, effective temperature is plotted 
versus the slope = d{6Aj''^)/d(xexc)- As can be seen, there is a 
clear linear relation between T^ff and the slope, with some minor 
spread of only 0.8 K due to a range in adopted surface gravities. 

Regarding the ionization equilibrium, we show in Fig. IC.2I 
the dependence between surface gravity and A^'~' (see eq. 3). A 
linear fit represents well this relation, with a scatter in log g of 
only 0.006 dex for models in a range of effective temperatures. 

In Fig. IC.3I we show the linear dependence between micro- 
turbulence velocity and the slope d{6A^'^^)/d(EWr)- For a given 
model, V, could be constrained to within 0.0004 km s ' . A range 
in Teff (5791 K < Teff < 5797 K) and log g (4.40 dex < log g < 
4.52 dex) imply in a scatter of only 0.009 km s"' in v,. 

Given the above dependences, the stellar parameters Tetj, 
log g and V, must be iteratively modified until the spectroscopic 
equilibrium conditions (equations 2-4) are satisfied simultane- 
ously. Since the degeneracy is relatively small, the final solution 
{rj/log g/v, = 5794 K/4.46 dex/l.OOkm s"') is very close to the 
independent solutions shown in Figs.|Cl]-|0](5794.5 K/4.458 
dex/1.006kms-'). 

In order to check how unique is the derived final solution, 
we have run over 200 models with different stellar parameters, 
with a very fine grid (steps of only 1 K in T^ff and 0.01 dex in 
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Fig. C.l. Tf.fl as a function of the slope = d(6Aj''^) / d(xexc) (see 
eq. 2), for a range of surface gravities (4.40 dex < log g < 4.52 
dex), corresponding to a scatter of only 0.8 K in Teff. The result 
of different model atmospheres are shown as open circles. The 
line represents a linear fit. 
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log g) near our best solution. We then evaluated how close to 
zero are the slope in T^ff (eq. 2) and the ionization equilibrium 
parameter A" (eq. 3). The following quantity is evaluated for 
each model. 



TG = (\slope /error] + |A'"Verror|)/2 



(C.l) 



The model showing the lowest TG value would be the best 
spectroscopic solution, which in our case is obtained for Tetj - 
5794 K, log g - 4.46 dex, and v, - 1.00 km s"'. A contour plot 
for the TG parameter is shown in Fig JC.4l Besides the best so- 
lution at TefF - 5794 K and log g = 4.46 dex, there are a few 
other nearby plausible solutions, with a mean value at Tetf - 
5794.3±0.5 K and log g = 4.462±0.012 dex, shown by a cross 
in Fig jC.4l Our grid samples a much larger coverage than that 
shown in Fig. IC.4I and we have verified that the best solution 
indeed represents a global minimum, i.e, there is no other solu- 
tion that can simultaneously satisfy the conditions of differential 
spectroscopic equilibrium. Thus, within the error bars our solu- 
tion is "unique". 

Appendix D: Helium abundance and the age and 
log g of HIP 56948 

Since our stellar parameters are very precise, actually the He 
abundance in HIP 56948 cannot be arbitrarily different from the 
solar He abundance. For example, an evolutionary track com- 
puted with a He abundance 5% higher than solar, would shift the 
Tetf by about -1-74 K at the same log g, i.e., a change 10 times 
larger than our error bar in Teff, leading thus to no plausible so- 
lutions. We are currently building an extensive grid of models 
with He as a free parameter, using the Dartmouth stella r evolu- 
tion code (IChabover et al.L 120011: iGuenther et al.L |1992|) . which 
is based upon the Yale stellar evolution code. 



In Fig. ID. II we show evolutionary tracks for M - 1.012 M© 
(solid lines), which is the best mass found for HIP 56948 using 
the Dartmouth tracks adopting the He solar abundance. These 
models were computed at three different helium abundances, so- 
lar and ± 1 % solar. We also show isochrones at 3 and 4 Gyr for 
different He abundances. The error bars in log g and Tgff put 
constraints on the He abundance, which can not be radically dif- 
ferent from solar Notice that the isochrones run parallel to each 
other and with only a minor shift for a change of ± 1 % in the He 
abundance, resulting thus in about the same central solution for 
age, independent of the adopted He content. Therefore thanks to 
our small error bars in stellar parameters we can put stringent 
constraints on the age of HIP 56948. Interestingly, there is a de- 
generacy between He and mass, although below our 2% error 
bar in mass. 

Another example where the adoption of a somewhat different 
He abundance did not affect much the stellar age is for the pair 
of solar analogs 16 Cyg A and B. For this pair, asteroseismology 
have recently constrained the ages of these stars, which are in 
excellent agreement with those derived from our isochrone tech- 
nique, despite the somewhat different adopted He abundances. 
Based o n three months of a lmost uninterrupted Kepler obser- 
vations, iMiiiaifcetaD (1201 2h obtained an age - 6.8±0.4 Gyr 
for their optimal models, in excellent agreement with an age = 
7.1±0.4 Gyr derived from our isochrone technique for the 16 
Cyg pair (Ramfre z et al.L 1201 ih . 

The effect of changing He by + 1 % (in Y) in our model atmo- 
spheres has als o a minor impact on the derived spectroscopic log 
g. As shown bv lStromgren et al1(ll982h . for solar type dwarfs the 
change in log g due to a change in the helium to hydrogen ratio 
(y = NhJNh) is: 



logg = logg' -Hlog 



(1+ 43; 0(1 +3') 



(l+4y)(l+y') 



(D.l) 
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Fig. C.3. Microturbulence velocity vs. slope d(SA^^^)/d(EW,-) 
(see eq. 4). The dotted, dashed and long dashed lines are for T^ff 
= 5791, 5794, 5797 K, respectively, and the spread shown for 
each line style is due to a range in log g (4.40 dex < log g < 
4.52 dex). This spread in T^ff and log g corresponds to a scatter 
of 0.009 km s ' in v,. The solid line represents a linear fit. 
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Fig. C.4. Contour plot of the parameter TG (eq. C.l), which 
evaluates how good is the differential spectroscopic equilibrium. 
The minimum is shown by a cross at T^ff - 5794. 3±0. 5 K and log 
g - 4.462±0.012 dex, which is in excellent agreement with our 
adopted solution. The contour levels increase in steps of ATG - 
0. 1 from the minimum. 



iLind etal](l2011Bh have shown that the relation above is also 
adequate for giant stars. For a change of + 1 % in Y in HIP 56948, 
the predicted change in log g is only -0.001 dex, which is well 
below our error bar in log g (0.02 dex). Thus, assuming that the 
He abundance of HIP 56948 is not radically different from solar, 
our derived spectroscopic log g value is essentially unaffected. 



Fig.D.l. Solar metallicity evolutionary tracks for 1.012 Mq 
(solid lines) at three different He abundances: solar and +1% 
solar. Isochrones at 3 and 4 Gyr are plotted with dashed lines. 
The position of HIP 56948 and error bars in T^ff and log g are 
also shown. 

Appendix E: v sin / and macroturbulence velocity 

We have determined v sin / from the differential line broadening 
(HIP 56948 - Sun). Naively we could be assuming an identical 
macroturbulence for both stars, but at a given luminosity class, 
macroturbulence seemsji_sniooth fu nction of temp erature (e.g., 
ISaar& OstenL 1 19971: iOravl l2005l IValenti & FischeR l2005i) . so 
neglecting this effect can lead to a slight overestimation of v sin 
i in HIP 56948 because it is hotter than the Sun and therefore the 
contribution of Vmacro to the line broadening in HIP 56948 should 
be slightly larger 

The tren d of macroturbulence velocity with reft described by 
iGravl (l2005l) for main sequence starsQ can be fitted by 

Vmacro = 13.499 - 0.00707 Teft + 9.2422 x lO^^r^^. (E.l) 

A similar correlation was advocated bv lValenti & Fischerl(l2005l) 
(after normalization to v®acro - 3.50 km s"', which is the value 
obtained for the Sun using Gray's relation): 

Vmacro = 3.50 + (T^eft - 5777)/650. (E.2) 

Finally, t he mean rela t ion (a ctive and non-active stars) ob- 
tained by 'Saar & Osten ( 1997), after transforming (B-V) to T^fi 
(IValenti & Fischer^ i2005i) and normalizing it to v®a„o = 3.50 km 
s"\ is: 

Vmacro = 3.50 + (Feft - 5777)/388. (E.3) 

The first two relations are valid for ~ 5000-6500 K, while the 
last relation is valid for ~ 5000-6100 K. On average, the above 
relations predict a differential (HIP 56948 - Sun) Avmacro = 
0.044 + 0.018 km s-i. 

For subgiants, giants an d luminous giants, our fits can be found in 
iHekker & Melende j l2007h 
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In order to determine v sin ; we selected 19 lines in the 602- 
682 run region, although essentially similar results are obtained 
(albeit with even so slightly larger errors) when 50 lines cover- 
ing the 446-682 nm region are used. First, we performed spectral 
synthesis of selected hues, in order to calibrate the relation be- 
tween line width (in A) and total broadening (in km s~'). Then, 
we estimated the total broadening using a much larger set of 
lines, and obtained vsin; after subtracting both the instrumental 
and the macroturbulence broadening. 

After taking into account the somewhat higher macroturbu- 
lence velocity of HIP 56948, we find v sin j/vsin = 1.006 ± 
0.014, or A V sin / = +0.013 ± 0.026 km s"' (or ±0.032 km s'^ 
including the error in macroturbulence), i.e., HIP 56948 seems 
to have about the same rotation velocity as the Sun, or rotating 
slightly faster, although it is unclear how much faster due to the 
uncertain sin i factor. 
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Table 1. McDonald radial velocity measurements for HIP 56948 



BJD 


Velocity 


Error 




(m s-') 


(m s-i) 


2454250.742247 


1.86 


3.51 


2454251.762237 


6.18 


6.37 


2454347.603043 


-0.26 


8.35 


2454822.022016 


-3.17 


4.48 


2455285.780573 


6.53 


6.00 


2455585.972936 


-5.06 


4.07 


2455643.770312 


-0.17 


4.67 


2455667.816687 


-16.11 


4.17 


2455990.845569 


10.21 


5.51 



Table 2. Keck radial velocity measurements for HIP 56948 



BJD Velocity Error 

(ms"') (m s"') 



2455610, 


,068981 


-2.77 


3.24 


2455611, 


,013099 


-6.06 


3.73 


2455611, 


,100610 


-9.49 


3.38 


2455611, 


,148864 


-10.18 


3.83 


2455766, 


,748023 


5.54 


2.20 


2455767, 


,742970 


4.39 


2.65 


2455767, 


,746153 


0.82 


3.78 


2455767, 


,749454 


3.74 


2.00 


2455935, 


,157734 


-0.35 


2.43 


2455935, 


,160951 


1.21 


1.89 


2455935, 


,163959 


0.67 


2.34 


2455936, 


,135676 


-4.52 


3.28 


2455936, 


,138532 


-5.07 


1.98 


2455936, 


,141412 


-8.27 


2.68 


2455937, 


,130793 


0.19 


2.85 


2455937, 


,133891 


-1.58 


2.46 


2455937, 


,137006 


-1.07 


2.54 


2455938, 


,111144 


3.23 


3.62 


2455938, 


,114027 


3.34 


2.50 


2455938, 


,116890 


-2.13 


4.05 


2455962, 


,089413 


1.70 


2.95 


2455962, 


,092257 


0.20 


2.62 


2455962, 


,095094 


0.85 


1.84 


2455963, 


,104952 


3.75 


3.12 


2455963, 


,107611 


5.20 


1.86 


2455963, 


,110271 


1.18 


2.34 


2455964, 


,037962 


1.92 


2.66 


2455965, 


,071634 


3.50 


1.98 


2455965, 


,075176 


4.35 


3.92 


2455965, 


,078440 


5.71 


2.31 



Table 4. Stellar parameters and Li abundance of HIP 56948 rel- 
ative to the Sun (HIP 56948 - Sun) 





Alogg 


A[Fe/H] 


Av, 


ALi (NLTE) 


Av sin / 


method 


reference 


(K) 


(dex) 


(dex) 


(km s-') 


(dex) 


(km s-') 






17±7 


+0.02+0.02 


+0.02+0.01 


+0.01+0.01 


0.23+0.05 


+0.01+0.03 


spectroscopy 


This work 


26±70 


-0.07+0.07 










IRFM, Hipparcos 


L. Casagrande (priv. communication) 


26+63 












IRFM 


Casagrande et al. (2011) 


24+25" 












IRFM 


Casagrande et al. (2010) 


17±5 


+0.01+0.01 


+0.02+0.01 


-0.01+0.01 


0.22 


+0.05 


spectroscopy 


Takeda & Tajitsu (2009) 


3±5 


-0.02+0.01 


+0.01+0.01 


-0.01+0.04 






spectroscopy 


Takeda & Taiitsu (2009) 


60±56 


+0.03+0.08 


+0.04+0.03 




0.22+0.07 


0.0+0.1 


spectroscopy 


Ramirez et al. (2009) 


5±36 


-0.04+0.05 


+0.01+0.02 


+0.01+0.06 


-0.02+0.13 


0.0+0.1 


spectroscopy 


Melendez & Ramirez (2007) 


-2±52 












photometry 


Masana et al. (2006) 



Notes. Error bar based only on photometric errors 
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Table 5. Atomic and molecular line list in the vicinity of the Li 
Unes 



Wavelength 


Species 




log gf 


A 




(eV) 


(dex) 


6706.5476 


CN 


3.13 


-1.359 


6706.5665 


CN 


2.19 


-1.650 


6706.657 


CN 


0.860 


-2.993 


6706.658 


CN 


0.614 


-3.622 


6706.728 


CN 


0.625 


-2.400 


6706.7329 


CN 


0.870 


-1.768 


6706.8440 


CN 


1.96 


-2.775 


6706.8626 


CN 


2.07 


-1.882 


6706.880 


Fe II 


5.956 


-4.103 


6707.00 


Si I 


5.954 


-2.56 


6707.172 


Fel 


5.538 


-2.810 


6707.2052 


CN 


1.97 


-1.222 


6707.272 


CN 


2.177 


-1.416 


6707.2823 


CN 


2.055 


-1.349 


6707.300 


C2 


0.933 


-1.717 


6707.3706 


CN 


3.05 


-0.522 


6707.433 


Pel 


4.608 


-2.25 


6707.460 


CN 


0.788 


-3.094 


6707.461 


CN 


0.542 


-3.730 


6707.4695 


CN 


1.88 


-1.581 


6707.473 


Smll 


0.933 


-1.91 


6707.548 


CN 


0.946 


-1.588 


6707.5947 


CN 


1.89 


-1.451 


6707.596 


CrI 


4.208 


-2.667 


6707.6453 


CN 


0.946 


-3.330 


6707.660 


C2 


0.926 


-1.743 


6707.7561 


'Li 


0.000 


-0.428 


6707.7682 


'Li 


0.000 


-0.206 


6707.809 


CN 


1.221 


-1.935 


6707.8475 


CN 


3.60 


-2.417 


6707.8992 


CN 


3.36 


-3.110 


6707.9066 


^Li 


0.000 


-1.509 


6707.9080 


^Li 


0.000 


-0.807 


6707.9187 


■'Li 


0.000 


-0.807 


6707.9196 


«Li 


0.000 


-0.479 


6707.9200 


'Li 


0.000 


-0.807 


6707.9230 


'Li 


0.000 


-0.178 


6707.9300 


CN 


1.98 


-1.651 


6707.970 


C2 


0.920 


-1.771 


6707.980 


CN 


2.372 


-3.527 


6708.023 


Si I 


6.00 


-2.80 


6708.0261 


CN 


1.98 


-2.031 


6708.0728 


"Li 


0.000 


-0.303 


6708.094 


VI 


1.218 


-2.922 


6708.099 


Ce II 


0.701 


-2.120 


6708.1470 


CN 


1.87 


-1.884 


6708.282 


Pel 


4.988 


-2.70 


6708.3146 


CN 


2.64 


-1.719 


6708.347 


Pel 


5.486 


-2.58 


6708.3700 


CN 


2.64 


-2.540 


6708.420 


CN 


0.768 


-3.358 


6708.534 


Pel 


5.558 


-2.936 


6708.5407 


CN 


2.50 


-1.876 


6708.577 


Pel 


5.446 


-2.684 
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Table 6. Stellar abundances [X/H] in LTE and NLTE, and errors 
due to uncertainties in the stellar parameters. 



Element 


LTE 


NLTE 


ATeff 


Alog^ 


Av, 


A[Fe/H] 


param" 


obs* 


total"^ 








+7K 


+0.02 dex 


+0.01 km 


+0.01 dex 










(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


C 


0.007 


0.007 


-0.004 


0.003 


0.000 


0.000 


0.005 


0.006 


0.008 





0.012 


0.011 


-0.005 


0.001 


-0.001 


0.002 


0.006 


0.008 


0.010 


Na 


0.014 


0.014 


0.004 


0.000 


0.000 


-0.001 


0.004 


0.006 


0.007 


Mg 


0.013 


0.012 


0.004 


-0.002 


-0.001 


-0.001 


0.005 


0.005 


0.007 


Al 


0.011 


0.010 


0.004 


0.000 


0.000 


-0.001 


0.004 


0.008 


0.009 


Si 


0.022 




0.002 


0.001 


-0.001 


0.001 


0.003 


0.003 


0.004 


S 


0.004 




-0.003 


0.003 


0.000 


0.001 


0.004 


0.005 


0.007 


K 


0.007 


0.007 


0.006 


-0.007 


-0.002 


0.001 


0.009 


0.009 
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Notes. Adding errors in stellar parameters Observational errors 
Total error (stellar parameters and observational) 
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Table 3. Adopted atomic data, equivalent widths, and differential NLTE corrections (HIP 56948 - Sun) 
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Table 3. Continued. 
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2.80 


7.9 


8.2 






5520.497 


21, 


.0 


1.8649 


0.55 


2.80 


6.8 


6.8 






5671.821 


21, 


.0 


1.4478 


0.55 


2.80 


15.6 


15.3 






4420.661 


21, 


I 


0.6184 


-2.273 


2.80 


16.2 


16.5 






5657.87 


21 


1 


1.507 


-0.30 


2.80 


69.4 


67.3 






5684.19 


21, 


\ 


1.507 


-0.95 


2.80 


39.5 


38.2 






6245.63 


21, 


\ 


1.507 


-1.030 


2.80 


34.2 


33.8 






6279.76 


21, 


\ 


1.500 


-1.2 


2.80 


30.8 


29.8 






6300.698 


21, 


1 


1.507 


-2.0 


2.80 


5.9 


5.6 






6320.843 


21, 


1 


1.500 


-1.85 


2.80 


9.9 


9.3 






6604.578 


21, 


1 


1.3569 


-1.15 


2.80 


37.1 


36.6 






4281.369 


22, 


.0 


0.8129 


-1.359 


0.502E-31 


24.6 


24.4 


-U, 


AA 1 
.UUI 


4465.802 


22, 


.0 


1.7393 


-0.163 


0.398E-31 


38.1 


37.9 


-U, 


.UUI 


4758.120 


22, 


.0 


2.2492 


0.425 


0.384E-31 


42.6 


42.9 


-U, 


AA 1 

.UUI 


4759.272 


22, 


.0 


2.2555 


0.514 


0.386E-31 


46.4 


46.3 


U, 


AAA 
.UUU 


5022.871 


22, 


.0 


0.8258 


-0.434 


0.358E-31 


70.1 


69.8 


u, 


AAO 
.UUZ 


5113.448 


22, 


.0 


1.443 


-0.783 


306E-31 


27.3 


27.6 


-U, 


AAO 

.UUz 


5219.700 


22, 


.0 


0.021 


-2.292 


0.208E-31 


29.3 


29.0 


A 


AAO 

.UUz 


5490.150 


22, 


.0 


1.460 


-0.933 


0.541E-31 


21.6 


21.4 


-U, 


AAQ 
.UU3 


5866.452 


22, 


.0 


1.066 


-0.840 


0.216E-31 


47.5 


47.9 


U 


AAO 

.UUz 


6126.217 


22, 


.0 


1.066 


-1.424 


0.206E-31 


22.1 


22.2 


A 

-U, 


AA1 
.UUI 


6258.104 


22, 


.0 


1.443 


-0.355 


0.481E-31 


51.9 


51.4 


A 

-U 


AA'2 


6261.101 


22, 


.0 


1.429 


-0.479 


468E-31 


49.5 


48.7 


-0, 


.003 


4583.408 


22, 


.1 


1.165 


-2.87 


2.80 


31.3 


31.4 


0, 


.000 


4636.33 


22, 


.1 


1.16 


-3.152 


2.80 


19.2 


19.1 


0, 


.000 


4657.212 


22, 


.1 


1.243 


-2.8 


2.80 


32.4 


32.2 


0, 


.000 


4865.611 


22, 


.1 


1.116 


-2.81 


2.80 


40.3 


39.3 


0, 


.000 


4911.193 


22, 


,1 


3.123 


-0.537 


2.80 


51.6 


51.8 


-0, 


.002 


5211.54 


22, 


.1 


2.59 


-1.49 


2.80 


34.2 


33.0 


-0, 


.001 


5381.015 


22, 


,1 


1.565 


-1.97 


2.80 


61.1 


59.9 


0, 


.000 


5418.767 


22, 


.1 


1.582 


-2.11 


2.80 


49.5 


48.8 


0, 


.001 


4594.119 


23, 


.0 


0.068 


-0.67 


0.216E-31 


55.6 


55.8 
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Table 3. Continued. 



Wavelength 


ion 


Xexc 






EW 


EW 


ANLTE 


(A) 




(eV) 






HIP 56948 


Sun 


(dex) 


4875.486 


23.0 


0.040 


-0.81 


0.198E-31 


45.9 


45.4 




5670.85 


23.0 


1.080 


-0.42 


0.358E-31 


19.5 


19.2 




5727.046 


23.0 


1.081 


-0.011 


0.435E-31 


39.4 


38.9 




6039.73 


23.0 


1.063 


-0.65 


0.398E-31 


13.8 


13.1 




6081.44 


23.0 


1.051 


-0.578 


0.389E-31 


14.7 


14.4 




6090 21 


23.0 


1.080 


-0.062 


398E-31 


33.8 


33.1 




6119.528 


23.0 


1.064 


-0.320 


0.389E-31 


23.4 


22.6 




6199.20 


23.0 


0.286 


-1.28 


0.196E-31 


14.6 


14.6 




6251.82 


23.0 


0.286 


-1.34 


0.196E-31 


16.2 


15.3 




6274.65 


23.0 


0.267 


-1.67 


0.194E-31 


9.2 


8.5 




6285.160 


23.0 


0.275 


-1.51 


0.194E-31 


11.2 


10.8 




4801.047 


24.0 


3.1216 


-0.130 


0.452E-31 


49.2 


49.4 


-0.003 


4936.335 


24.0 


3.1128 


-0.25 


0.432E-31 


44.3 


44.3 


-0.003 


5238.964 


24.0 


2.709 


-1.27 


0.519E-31 


16.5 


16.3 


-0.002 


5247.566 


24.0 


0.960 


-1.59 


0.392E-31 


81.3 


81.1 


-0.003 


5272.007 


24.0 


3.449 


-0.42 


315E-30 


24.6 


24.0 


-0.003 


5287.20 


24.0 


3.438 


-0.87 


0.309E-30 


11.4 


11.5 


-0.003 


5783.09 


24.0 


3.323 


-0.43 


0.802E-30 


31.8 


31.2 


-0.002 


4588.199 


24.1 


4.071 


-0.594 


2.80 


70.8 


69.7 


-0.003 


4592.049 


24.1 


4.073 


-1.252 


2.80 


47.5 


47.7 


-0.002 


5237.328 


24.1 


4.073 


-1.087 


2.80 


54.0 


54.0 


0.000 


5246.767 


24.1 


3.714 


-2.436 


2.80 


16.1 


15.5 


0.001 


5502.067 


24.1 


4.168 


-2.049 


2.80 


19.3 


18.7 


0.000 


4082 939 


25.0 


2.1782 


-0.354 


255E-31 


90.0 


89.8 


-0.007 


4739.10 


25.0 


2.9408 


-0.490 


0.352E-31 


61.5 


60.3 


-0.002 


5004 891 


25.0 


2.9197 


-1.63 


0.314E-31 


14.0 


14.1 


-0.004 


6013.49 


25.0 


3.073 


-0.251 


2.80 


87.8 


86.5 


-0.005 


6016.64 


25.0 


3.073 


-0.084 


2.80 


97.0 


96.1 


-0.006 


6021.79 


25.0 


3.076 


+0.034 


2.80 


89.5 


89.3 


-0.005 


5212.691 


27.0 


3.5144 


-0.11 


0.339E-30 


21.4 


20.6 


-0.002 


5247.911 


27.0 


1.785 


-2.08 


0.327E-31 


18.1 


17.7 


-0.001 


5483.352 


27.0 


1.7104 


-1.49 


0.289E-31 


51.1 


51.2 


-0.002 


5530.774 


27.0 


1.710 


-2.23 


0.226E-31 


18.9 


19.1 


-0.002 


5647.23 


27.0 


2.280 


-1.56 


0.414E-31 


15.0 


14.7 


-0.002 


6189.00 


27.0 


1.710 


-2.46 


0.206E-31 


11.6 


11.3 


-0.002 


6454.995 


27.0 


3.6320 


-0.25 


0.378E-30 


14.7 


14.1 


-0.002 


5589.358 


28.0 


3.898 


-1.14 


0.398E-30 


27.6 


27.0 




5643 078 


28.0 


4.164 


-1.25 


379E-30 


16.1 


15.7 




6086.282 


28.0 


4.266 


-0.51 


0.406E-30 


45.7 


44.9 




6108.116 


28.0 


1.676 


-2.44 


0.248E-31 


64.3 


63.3 




6130.135 


28.0 


4.266 


-0.96 


0.391E-30 


24.4 


23.6 




6204 604 


28.0 


4.088 


-1.14 


277E-30 


23.1 


22.4 




6223.984 


28.0 


4.105 


-0.98 


0.393E-30 


28.7 


28.3 




6378.25 


28.0 


4.1535 


-0.90 


0.391E-30 


32.5 


31.8 




6161.112 


28.0 


1.826 


-2.17 


2.80 


79.9 


79.2 




6772.315 


28.0 


3.657 


-0.99 


0.356E-30 


50.6 


50.0 




7727.624 


28.0 


3.678 


-0.4 


0.343E-30 


91.4 


90.7 




7788.930 


28.0 


1.950 


-2.0 


0.218E-31 


92.7 


92.1 




7797.586 


28.0 


3.89 


-0.34 


2.80 


80.0 


78.6 




5105.541 


29.0 


1.39 


-1.516 


2.80 


91.5 


91.5 




5218.197 


29.0 


3.816 


0.476 


2.80 


52.2 


51.2 




5220.066 


29.0 


3.816 


-0.448 


2.80 


17.3 


17.3 




7933.13 


29.0 


3.79 


-0.368 


2.80 


30.5 


30.8 




4722.159 


30.0 


4.03 


-0.38 


2.80 


71.9 


71.2 


-0.001 


4810 534 


30.0 


4.08 


-0.16 


2.80 


73.8 


73.3 


-0.001 


6362.35 


30.0 


5.79 


0.14 


2.80 


21.6 


21.1 


0.000 


4854.867 


39.1 


0.9923 


-0.38 


2.80 


48.3 


48.4 




4883.685 


39.1 


1.0841 


0.07 


2.80 


57.6 


57.1 




4900.110 


39.1 


1.0326 


-0.09 


2.80 


55.5 


55.3 




5087.420 


39.1 


1.0841 


-0.17 


2.80 


49.1 


48.4 




5200.413 


39.1 


0.9923 


-0.57 


2.80 


39.0 


38.6 




4050.320 


40.1 


0.713 


-1.06 


2.80 


23.8 


23.4 


0.001 


4208.980 


40.1 


0.713 


-0.51 


2.80 


43.2 


42.0 


0.000 


4442.992 


40.1 


1.486 


-0.42 


2.80 


25.7 


24.7 


0.000 


5853.67 


56.1 


0.604 


-0.91 


0.53E-31 


64.6 


63.7 


-0.002 


6141.71 


56.1 


0.704 


-0.08 


0.53E-31 


116.2 


115.7 


0.000 
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Table 3. Continued. 



Wavelength 


ion 


Xexc 


log^/ 


Ce 


EW 


EW 


ANLTE 


(A) 




(eV) 






HIP 56948 


Sun 


(dex) 


6496.90 


56.1 


0.604 


-0.38 


0.53E-31 


99.5 


99.2 


0.000 



